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Chapter 1 : Introduction 

The past few decades have seen helicopters evolve into 
truly multi-purpose machines. In the cities, radio stations 
use them to provide- motorists with valuable • traffic 
information, and in some of the largest they have become just 
one more form of commuting. In the battlefield, helicopters 
have proved themselves effective tools in rescue operations, 
as well as for the transport of troops and supplies to 
otherwise inaccessible regions. In either role, civilian or 
military, they are expected to do their job as quietly as 
possible. 

Everyone has heard the unpleasant "slapping” noise from 
an approaching helicopter. When present, and it often is, in 
hover, it usually dominates all other sources of helicopter 

noise - including engine noise. Investigators generally 

( 1 ) 

agree'- ‘ that for moderate subsonic tip Mach numbers a great 
part of each slap signal has its origin in the impulsive load 
a blade feels upon passage over a vortex trailing from 
another. For higher tip speeds, however, Boxwell, et, al. 
have demonstrated experimentally that blade slap can be 
mostly due another aerodynamic effect: the presence of 

tip-region shocks on the advancing side of the rotor. Here 
we develop aeroacoustic models to predict blade slap due to 
blade-vortex interaction alone. 

Normally, heavier helicopters operate ,f,aster turning 
rotors for greater values of needed lift (instead, the number 
of blades can sometimes be increased). The analysis to 
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determine the unsteady blade force and associated acoustic 
field for a blade-vortex interaction at high subsonic 
tip-Mach numbers is fundamentally different from that for low 
tip-Mach numbers; compressibility effects then . become 
important in the calculation of a section of the signal’s 
spectrum which may contain a substantial fraction of the 

signal's acoustic energy. Also, compressibility changes the 
qualitative character of acoustic radiation: pulses from 

blade surface dipoles firing at different values of retarted 
time t-r/Cg may nevertheless reach the far-field point r 
simultaneously at time t; for high subsonic M, these, 
differences in firing time A(t-r/Cg)'s in general are not 

small in comparison to the period of oscillation of 
high-level spectral components of the total slap signature. 
In the present study we develop models to predict 
theoretically this important part of the signal's spectrum, 
the part which causes our blade to behave largely as a 

noncompact acoustic source. We shall see that in the 

boundary-value problems we need to solve compressibility and 
high-frequency effects appear lumped in two parameters, both 
proportional to reduced frequency divided by 

In the linear models developed here, we split the 
aeroacoustic problem of blade-vortex interaction into a sum 
of simpler gust problems; each gust represents a component of 
the spectrum of the vortex-induced upwash on the plane of 

i • 

flight. It is important to point out that application of our 
models for blade-gust interaction is not limited to 
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predictions of noise due to blade-vortex interaction, but 
that they could also be used, for example, to calculate 
unsteady forces acting on compressor blades in the presence 
of nonuniform incident flows, or to predict the acoustic 
field surrounding a. thin wing cutting through confined 
regions of turbulence. 

In earlier theoretical studies of helicopter-blade 
slap,^^’"*^ investigators have applied incompressible flow 
models to calculate unsteady sectional loads acting on a 
blade passing over a vortex. They used this computed lift 
force in three-dimensional acoustic models to represent the 
strength of simple point dipoles along a chosen portion of 
the blade length. Application of such chordwise-compact 
aeroacoustic models of blade-vortex interaction to predict 
noise for practical cases of interest was justified whenever 
interaction parameters met one of the following two 
conditions: (1) Large vertical blade-vortex separation to 
chord ratio for moderate Mach numbers; (2) very low Mach 
numbers for arbitrary blade-vortex separations. 

Actual helicopter blades often operate such that neither 
of these two compactness criteria is satisfied, however, and 
noncompact aeroacoustic theories such as those developed here 
become necessary. The natural uncoupling that exists between 
an acoustic field and its aerodynamic source (in 
two-dimensions, the loading along the chord) disappears when 
the latter is acoustically noncompact; there is no 
incompressible flow region surrounding the airfoil then, and 
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the acoustic field begins at the airfoil surface. Instead of 
radiation, from a simple dipole, we expect a more complicated 
acoustic field due to a distribution of dipoles over the 
chord with large differences in phase. Kaji has discussed in 
detail noncompactness effects in ref. 5. 

Only a few of many important three-dimensional problems 
in unsteady aerodynamics have been solved in closed form for 
subsonic flow; for example, until very recently an unsteady 
lifting-line theory was not available (ref. 6). Another 
three-dimensional, but chordwise- compact, model is that due 

(7 3 ^ 

to Chu' ’ , who developed a numerical lifting-surface theory 
for a semi-infinite blade passing through a gust, and then 
applied it to predict noise due to blade-vortex interaction. 
Many past efforts in modelling the strength of acoustic 
sources due to blade-vortex interaction, however, have been 
based on two-dimensional aerodynamic theories. These models 
assumed that at most spanwise positions on a blade 
three-dimensional tip effects can be neglected; they become 
more and more valid as the vertical blade-vortex separation 
is reduced and the loading at each section along the length 
of the blade is dominated two-diraensionally by local 

chordwise unsteadiness. 

Cl 3 ) 

Landahl has shown that the surface pressure 
distribution for a thin two-dimensional airfoil passing 
subsonically through a gust of arbitrary wavelength may be 
determined by means of an iterative scheme applied in the 
aerodynamic reference frame (fixed on the airfoil). In the 


t 



first step the chord is allowed to extend infinitely in the 
downstream direction from the leading edge; this eliminates 
the Kutta condition and the wake from the problem. The 
solution to this leading-edge problem satisfies both the 
upstream boundary condition <5 sO and the flow tangency 
boundary condition on the airfoil. A trailing-edge problem 
is then solved in order to correct the leading-edge solution 
in the wake. The sum of these two results satisfies flow 
tangency on the airfoil and pressure continuity across the 
wake; however it violates the condition of $ =0 upstream. 
Further iterations may be carried out to provide a more 
accurate solution. 

Landahl^^^ has shown that the series which results from 
aplication of his iteration scheme converges uniformly; but 
that due to the strong communication that exists between 
leading and trailing edges when the gust wavelength is much 
longer than the airfoil chord, many terras would be needed to 
obtain a good approximation to a solution which is easily 
available by. simpler methods because the source region is 
then acoustically compact. Landahl's iteration scheme 
becomes useful when the gust wavelength is much shorter than 
the airfoil chord. Then the leading edge and the trailing 
edge become largely independent of each other, and only a few 
terms in his series are needed to approximate a solution 
which is, difficult to obtain by any method known because the 

i • 

source region is acoustically noncorapact. In the 
short-wavelength limit the solution is given by the first 



term in the series; the relatively simple result from the 

leading-edge problem. 

( 12 ) 

Adamczyk has also given an infinite-series solution 
for a two-dimensional airfoil passing through a. gust of 
arbitrary wavelength. His approach was quite different to 
that of Landahl’s and the terms in his series of Mathieu 
functions cannot be interpreted physically as edge effects. 
He presented numerical results for some values of the 
acoustic directivity in the aerodynamic reference frame 
showing enhancement of the acoustic, signal in the downstream 
direction. 

fl3) 

Adamczyk' has also determined the response of an 
infinite-span swept wing to an oblique gust convected by a 
subsonic freestream. His expression for the pressure 
distribution on the wing was the sum of the first two terms 
of the series from the iteration scheme described above. 
Although he presented numerically calculated loads on an 
airfoil passing over a potential vortex in this work, he did 
not report any acoustic results, not even for the gust 
problem. 

(14) 

Amiet gave numerical results for the power-spectral 
density of the three-dimensional acoustic far field in the 
aerodynamic reference frame due the interaction of 
small-scale turbulence with a high-speed thin wing of finite 


span. He replaced the 

surface 

of 

the wing 

with 

a 

distribution of point 

dipoles 

of 

f • 

strengths 

given 

by 

Adamczyk' s two-dimensional response 

function 

for 

a 



short-wavelength gust, which he then used in a spanwise 
superposition to model finite-span effects. His 
three-dimensional acoustic model for numerical calculations 
is valid in- the high-frequency limit, or alternatively, in 
the limit of infinite aspect ratio. Amiet has estimated 
the power-spectral density of the acoustic far field due to a 
fast climbing rotor cutting through small-scale turbulence. 
The small-scale character of the turbulence allowed him to 
consider spanwise blade elements as moving in a locally 
rectilinear fashion. The results were given in a coordinate 
system fixed to the rotor hub. 

(IS) (13) 

Amiet has also generalized Adamczyk's original 
two-dimensional high-frequency theory by considering a gust 
not convected by the freestream. In addition to the gust 
problem, in this work he also investigated plunging motion of 
the airfoil; no acoustic analysis was performed, however. 

In chapter 2 we develop a two-dimensional aeroacoustic 
model for a blade of infinite span passing through an oblique 
gust of short wavelength. We apply Landahl’s concept of edge 
separation and obtain an expression for the pressure field 
due to the airfoil-gust interaction. The result is 
approximate because it is the sum of the first two terms in 
the iteration scheme of edge separation; the first 
leading-edge solution; and the first trailing-edge solution. 
Both of these two problems are ' solved in the aerodynamic 
reference frame by the VJiener-Hopf technique. The total 
solution contains the aerodynamics at the airfoil surface, 



that is, the surface pressure distribution, and the acoustics 
of the interaction in a rather compact analytic expression. 
Letting the field point in the solution come down on the 
airfoil, we recover Adamczyk's two-dimensional aerodynamic 
result for the case of zero sweep. Letting the field point 
go to infinity, we obtain the two-dimensional acoustic field 
in the airfoil reference frame due to the interaction of an 
oblique short-wavelength gust with a subsonic airfoil. This 
section of our presentation is a review of part of the work 
reported in ref, 9. 

In chapter 3 we perform a spanwise superposition of the 
two-dimensional solution from chapter 2 to approximate the 
three-dimensional pressure field for a rotating rectangular 
wing with subsonic tip-Mach number passing through a gust; 
this last result is given in the acoustic, or ground, 
reference frame. We calculate the three-dimensional acoustic 
far field by allowing the observer’s coordinates to go to 
infinity, and obtain a closed-form expression for the 
acoustic directivity pattern which does not require numerical 
computation . 

In ref. 9 the author performed a spanwise superposition 
of the two-dimensional result of chapter 2 which is different 
from the one presented here. The earlier one in ref. 9 is 
actually valid only in the infinite frequency limit, whereas 
the one we develop here in chapter 3 is valid for arbitrary 
high frequencies, not just the infinite limit. As we will 
see, the analysis in ref. 9 hid some interesting physical 
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features of the solution for cases in which the disturbance 
through the fluid has a supersonic trace speed. Also, the 
three-dimensional model developed in ref. 9 used a "box" 
spanwise superpositipn in which inboard blade sections 
reached the same maximum value of loading as sections near 
the blade tip. Here we develop a much better model for the 
unsteady loading on a rotating blade, which in reality sees a 
freestream linearly increasing from hub to tip: We model the 
three-dimensional unsteady blade load by the two-dimensional 
value at the tip, and say that it decreases from this value 
linearly to zero at the blade hub; the Mach number in the 
two-dimensional solution is used as the tip-Mach number in 
the three-dimensional solution. 

The actual loading on a helicopter blade, however, 
vanishes at the tip as /L-y, a functional behavior with 
infinite slope at y=L, where we therefore expect to find 
generation of strong trailing vorticity. It follows that our 
aerodynamic superposition model is not consistent with 
reality near the tip: it overestimates the strength of 
acoustic dipoles there, and also gives a wrong representation 
for their phase. As we discussed earlier, however, at high 
frequencies the spanwise extent of aerodynamic influence of 
the tip is reduced; blade sections react then mainly to 
strong local shed vorticity. Acoustic fields predicted by 
the model in chapter 3 s"hould be good approximations to those 
we would obtain using a model with correct magnitude and 
phase of tip dipoles. 
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( 17 ) 

Lima- has applied the three-dimensional acoustic 
theory developed in ref. 9 in a study of effects of 

(la) 

tip-vortex structure on radiated noise. He used Betz's 
theory to construct vortices from assumed forms of steady 
loading on a helicopter blade, and then computed numerically 
the spectrum of the upwash which each vortex investigated 
induced on the plane of flight. 

In the second section of chapter 3 we apply the 

leading-edge part of the acoustic theory developed in the 

/ 

first section to calculate the noise signature for the 

interaction with a potential vortex. The result is given in. 

closed form. It may be used also to -predict noise due to the 

interaction with a real vortex, which has a viscous core. 

(“) 

This follows from Widnall's suggestion that the aerodynamic 
effect of a real vortex may be modelled by a potential vortex 

at a somewhat larger effective distance below the flight 

plane. 

We also should point out here, as we did in ref. 9, that 
results equivalent to those obtained by our approach could be 
reached by an alternate method: By using the airfoil surface 

pressure distribution already available as the strength 
of point dipoles in a Green's function integral. With the 
present technique, however, the whole pressure field is 
obtained at once, so that the solution is both more elegant 

and compact. Also, the results here in terms of contour 

/ • 

integrals make the far-field calculations straightforward and 
unambiguous. For example, the phases of- certain quantities 
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are always obvious by the present complex-variables method. 
Finally, our solutions at a general field point may be used 
immediately to calculate other off-fairfoil flow properties of 
interest besides the acoustic field - the shed vorti.city in 
the wake, for instance. 

We have said earlier that for very short wavelengths the 
blade loading is well approximated by the two-dimensional 
result only at a distance away from the tip. It turns out, 
however, that in single-rotor helicopters blade-vortex 
interactions frequently occur in which only outboard sections 
of each blade are actually expected to radiate noise. The 
purpose of chapter 4 is to develop a three-dimensional 
aeroacoustic model for a blade tip passing through a gust. 
The objective is to determine . the actual magnitude and 
spanwise phase variation of tip dipoles. Such a tip model 
should replace chapter 3*s whenever conditions of 
noncompactness prevail, but frequencies are not high enough 
to ignore three-dimensional tip effects. 

Due to the complexity of the problem in general, we 
consider here the case of a square tip with side edge. Also, 
since earlier studies have demonstrated that in very 
noncompact situations the trailing edge can be neglected, we 
simplify the model further: We remove the trailing edge and 
determine the unsteady pressure field for a quarter-inf inite 
plate with side edge passing through a gust at high subsonic 


speed. 


Miles 




solved exactly a problem similar to the one we 
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consider here, but for supersonic flow. For the subsonic 

flow case no simple procedure yielding an exact closed-form 

result is available. In chapter 4 we obtain an approximate 

solution analytically for the pressure field everywhere in 

( 20 ) 

space; we apply a procedure similar to that Carrier Uas 
used to solve a three-dimensional diffusion problem. By 
allowing the field point down on the plate's surface we are 
able, with some work, to obtain a relatively simple 

expression for the tip loading valid for high-frequencies and 
for large values of the interaction angle A. In fact, this 
load expression turns out to be simple enough to enable us, 
in chapter 5, to reintroduce the trailing edge into the 
analysis and so arrive at an aerodynamic model for a 
rectangular-tip blade of semi-infinite span, with /y 

side-edge loading behavior, passing through an oblique gust. 
We then proceed to calculate the associated acoustic field in 
closed form. 

We end this introductory chapter with a very brief 

qualitative discussion on how we should expect a typical 

( 21 ) 

blade-slap signal to look. General acoustic theory tells 

us that a blade-slap pulse is proportional to the first time 
derivative of the unsteady blade force induced by the 
interaction. Fig. 1.1a shows a rotor turning at angular speed 

; the position of rolled-up tip vortices is indicated. 
Consider the particular blade-vortex interaction for which we 
have drawn the symbol A , for the interaction angle; The 
position of the vortex relative to the blade at that instant 



corresponds to time 1 in fig. 1.1b. Point 1 in fig.l.lc 
indicates the value of unsteady blade force F(t) at time 1; 
because its positive part nearly cancels its negative part 
there, the total load magnitude is relatively small. . At time 
3, however, F(t) reaches its largest negative value, and its 
time derivative dF/dt quickly shoots up as the vortex 
suddenly crosses outboard of the tip. At time 4 the net 
force is small again. 
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Chapter 2: Two-dimensional aeroacoustic theory for 
noncompact blade-gust interaction 

In this chapter we apply the concept of separation of 
leading- and trailing-edge effects at high reduced 
frequencies to obtain an approximate solution for the 
two-dimensional pressure field for an infinite-span wing 
passing through an oblique gust. The result is the sum of 
the solutions of the first leading- and ' trailing-edge 
problems. Before using it to calculate the two-dimensional 
acoustic field, we bring the field point down to the airfoil 
surface and recover the expression in ref. 13 for the 
pressure distribution due to the interaction with the gust; 
this suggests that our result for the pressure at a general 
field point in space must be correct. We then proceed to 
calculate the acoustic field by letting the field point in 
the solution go to infinity. The resulting expression for 
the directivity of sound is our contribution to the 
noncompact two-dimensional theory of airfoil-gust 
interaction; it first appeared in ref. 9. 

( 15 ) 

Following Amiet's formulation, we have in the 
aerodynamic reference frame a rigid flat-plate airfoil on the 
x-y plane with its chord extending from x=0 to x=2. It 
interacts at an angle A with a short-wavelength sinusoidal 
gust convected by the freestream (fig. 2.1). The linearized 
equation for the perturbation velocity potential $ is 


XX 


+ $yy + <I> 


ZZ 




Dt" 


0 


( 2 . 1 ) ; 
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with boundary conditions 


$(x<o,y,o,t) = 0 

(2.2a), 

iuit 

(2.2b), 

$^(o<X<2,y,o,t)= -bw(x,y) e 

^ (x>2,y,o,t) = 0 

(2.2c); 


where the spatial variables have been nondimensionalized with 
respect to the airfoil semichord b. D/Dt denotes 
9/3 t+(U/b) 3/3 X, the linearized substantial derivative. For a, 
sinusoidal gust of small amplitude w^ , the gust downwash 
w(x,y) becomes w^exp[i(o)t-kj^x-kyy)] ; for one convected by the 
freestream, the nondiraensionalize.d gust wavenumber in the x 
direction k^^ is then cjb/U, the reduced frequency, and the 
spanwise wavenumber ky=(ub/U) tanA . Boundary condition (2.2b) 
is a statement of flow tangency on the airfoil surface. From 
the linear relation • between potential and pressure, i.e., 
P=- p^(D$/Dt ) , boundary condition (2.2c) requires that there 
be no pressure discontinuities at the trailing edge and 
across the wake. 

Since the plate is infinite in the y direction, the 
potential $ will have the y dependence of the input gust, a 
traveling sinusoid; we therefore assume a harmonic behavior 

iojt-ikyy /,'• 

$(x ,y,z,t)=<{'(x,z) e ^ 

according to (2.2b). We also make the change of dependent 
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variable 

(j)(;^z)= <ji*(x, 2 ) e 

and compress the z coordinate by the change /I -M^ z=Z, The 
object of these last two transformations is to reduce the 
relatively complicated convected wave equation to the more 
familiar Helmholtz equation. 

The boundary-value problem for cji* becomes 

“0 (2.3), 

XX 22 


with 

I 


({,* (x<0,0) = 0 


(2.4a) , 


(o<X<2,0) = - bwQ a 



(2.4b) , 


(2.4c) ; 


where 


kj 

1 -^ 


^ - 




tan^A 


-1 

_k| 

-tan^A 


s i n^A 

" wF 

1-M^ 



- 


■* ^ 


J(2.4d), 


and is one of the two wavenumber-related similarity 
parameters in the problem; the other is kx/(1-M^). The 
parameter is negative for M/sinA<1 ( tanA>M/ /i - m-^ } . This 


case 


corresponds physically to a subsonic speed of 
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propagation of the disturbance due to the passage of the 
infinite-span wing^ through the stationary gust. Since a 
subsonically traveling disturbance of infinite extent 
radiates no sound, we expect no sound radiation for those 
Mach number and interaction angle combinations for which 
is negative. 

Leading-edge problem 

(0 

The boundary- value problem for (|>* , the velocity 

potential field due to the interaction of the gust with a 
semi-infinite chord airfoil extending downstream from the 


leading edge, is 

^*xx ^*ZZ 

+ =* 0 

(2.5); 

with 

4’^ (x<0,0) 

= 0 

( 2.6a) , 


41^ (x>0,0) 

a -bw-i exp [ “ikyX 

_] 


•*z 


(2.6b) . 

The problem 

may be readily 

solved by the 

Wiener-Hopf 

technique. We define the transform pair 




f • 

(2.7a) , 


? (x,7) = r 

dx^ 

rZlT 

(2.7b) ; 



where 


?(X,2)=$^a.Z) •+ 

that is, 

dx 4>*(x,2) 

*' (2.8a), 




Jo /2? 


(2.8b) . 


The contour C and the regions of analyticity 0, 0 in the 
complex X plane are yet to be determined by the physical 
requirement that acoustic waves propagate away from the 
airfoil, the Sommerfeld radiation condition. Transforming 
(2.5) and solving the resulting equation, we obtain that for 
this lifting problem 


5 U,Z±)= ± ? 


Z /X^ 


( 2 . 9 ) ; 


where the branch of is chosen so that its argument 
vanishes as lx| goes to infinity along C. We initially take 
p as a positive constant corresponding to M/sinA>l. The 
results may be analytically continued later to include the 
negative p^ case. 

Since an exp(iut) dependence was assumed, outward 
propagation requires that k^, and thus also ky and p , have a 
small negative imaginary part. The regions of 0 and 0 



analyticity in the complex X plane for (2.8a, b) are then 
given by Im(X)>Im(u), Im(X)<Ira(-y), respectively, as shown in 
Fig. 2.2. The contour C runs below the real axis for 
Re(X)<Re(-y) and above the real axis and the point 
for Re(X)>Re( y) . 


From (2.6a, b) and (2.8a, b) the following Wiener-Hopf 
equation is obtained (0 = 0) 



the left side of (2.10) being a Q function and the right a0. 
They are, therefore, analytic continuations of each other and 
so they are both at least entire functions of X, If 4 1 Z 
are assumed to be well behaved a X= ®, (the Riemann-Lebesgue 
lemma) Liouville's theorem then states that both sides of 
(2.10) are independently equal to zero. The right side of 
the equality says that 


? (X,o‘^)-l^o 1 


/ 2ir' 






( 2 . 11 ); 


and so, from (2.9) and (2.7a) 



(x.Z) = }bwo 
2Tr/17i^ 





-iXx - |Z I 
’ dX e 1 

I (X- ) 

^ l-W 


( 2 . 12 ) . 
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Arguments similar to those made for , the time and y 
dependence of the velocity potential also apply to the 
pressure. Making the same changes of dependent and 
independent variables^ made before for the potential, we have 


P(x,y,2,t) = p(x,z) 
p(x,z) = p*U.z) sxp [ ' ] 


(2.13a), 
(2.13b) . 


3 


I 


f 




I 


! 


From the relationship between pressure and velocity 
potential, the pressure p*^^\x,Z) corresponding to (j)^^(x,Z) is 
obtained ; 


. 0 ) 


(x,Z)= - OnWn ^ . V I Z 
2ir/T 




dX e 
>^X-u 


i Xx- j Z [ /X^-p^ 


(2.14) . 




For x>0, Z=0+, the contour C may be deformed to C in 
fig. 2.3 so that 


(x>0,oi) a + PqWq u 


/I 


1-M2 


+ U 



(2.15) , 


in agreement with the aerodynamic result in ref. 13. 
Trailing-edge problem 

The solution to the trailing-edge problem is a pressure 

pC.^)such that pj^^)+p|^^) is zero at the trailing edge and in the 

( ■ 

wake. From the linear relation between potential and 
pressure, we note that i^(^)(x,Z) also satisfies Helmholtz's 


1 


I 



equation. We also require that .3 p; V3Z vanish upstream of 

the trailing edge and that it be continuous for all x. With 

x=x'+2, the trailing-edge correction boundary-value problem 
(2) 

for p* 'in the x'-Z coordinate system is 


with 


(21 . ( 2 ) .2 ( 2 ) 

p*x'x' U P* 


0 


(2.16) ; 


3Z 


(x‘+2,0) = 0 for x'<0 


(2.17a), 


(x’ + 2 , 0 ^) » - P^^^ (x«+ 2 ,oi) fopxSo (2.17b). 

This problem may also be solved by the V/iener-Hopf 
technique. We define the transform pair 


(-V.z) 


P (X,Z) = 


^ P(X.Z) 

C 


^ ° p1^’ (?<' .2) 

i'Tn 


Where p= p^ p^ that is, 

P0(x',Z) j' Pp^CpZ) 


(2.18a) , 


(2.18b) ; 


(2".i9a), 


— <A 


?@CX.Z) = Pp’ (C.2) 


(2.19b). 


Transforming (2.16), solving the resulting equation and 
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( 2 ) 

invoking the continuity of 3 p* f for all x', we obtain 


pU.Z-) = ± P ® 




--y‘ 


( 2 . 20 ) ; 


where P<\ ,0 +) sPqCx ,O+) + P 0 ( x,0+) , the arbitrary function of X 
to be determined by the Wiener-Hopf technique. 

From boundary condition (2.17a), we see that 
3 P 0 ^X ,0)/3Z=0 . From boundary condition (2.17b), and 
.recalling the result in ( 2 . 15 ), ? 0 ( may be calculated; 


P (X. 0 -) U 

^ ]/2 


( 2 . 21 ) 


dt e 


it 


•a Fresnel integral; EC2(X-vi)] therefore 


where E(a)= 

Jo ’^2TTt 

has a branch point at X=p. 

The Wiener-Hopf equation in this case is ( 0 = 0 ) 

in/k 


p0(X,O+)+K' 


:* /X+U -2IX ( e' 


- E [2(X-y)] 


} 


■iTr/4 


tt /2 




. 1 — (X,0+) - K* (2-22); 


where 


K* = (1qWq_ (j g 


•\Tt/k 


■M 2 


+ y 


/ 1 -M 2 


a constant, and K f 0 denotes that additive part of /X+y 
P^X>0 + J which is analytic in the upper half x plane. The 
standard procedure of splitting an arbitrary function into a 
sum of two, one 0 and the other 0 ‘, is treated in detail in 


1 
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ref. 22. Making the usual assumptions regarding the good 
behavior of I^(X,0 + ) and 3l^(X,0)/3Z, we have that 
(Liouville's theorem) both sides of (2.22) vanish 
independently. This gives I^(^,0+). Adding (2.21) to it, we 
obtain P(A,0+), and substituting into (2.20) and by means of 
(2.18a), we finally have that 


pi^^(x.Z) ’-iPs^aJL- llM- 


^g-U(x-2)- \L I 
C 


d£ ^-2ic 

ftXVrrvre (2.23); 


where the contour C is indented above the real axis at X= C. 
Total Solution 


Adding (2.14) and (2.23) and applying (2.13b), we obtain 
the two-term approximation to the two-dimensional pressure 
field for the passage of an infinite-span wing through an 
oblique gust 


p(x,2) = P^'>(x.2) + (x,2) (2.24). 


\ 


The pressure distribution on the airfoil, with its chord 
extending from x=0 to x=2, is calculated from (2.24) which 
now says that 

p (0<x<2,0+) = (x>0,0'^) + p^^^ (^2,0'') (2.25). 

(i) 

The part p (x>0,0+) has already been determined and is 

( 2 ) 

given by (2.15) and (2.13b). To obtain p (x<2,0+), we start 
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with (2.23) and observe that for Z=0+, x<2, the integral 


• 00 ■ 

I(y,x,Z^) = f dX 

JQ Jh 


(2.26) 


may be deformed to 


^ -iX(x-2) r /OTe'^-^ 

l(u,x,0+) = d^e- (2.27); 

A0.4, . ^ 'S f- 


where the contour C**is shown in fig. 2.4. 

Interchanging orders of integration we have 


l(y,x,0+) = - d? d±e'' 


-iX(x-2) 


(2.28). 


The term 1/(X-5) attains its. maximum magnitude of -1/2U 
when ?= u and X=-u. and goes to zero as or as X-> -» along 

iHe 

C , or both. We can expect then that the largest 
contribution to l(u»x,0) comes from values of C near the 
lower limit. We therefore make the following approximation: 


I (P.x.O*) = - ^ e-'*' 


oo 

dCe'2''^f -I(X-u)(x-2) 

1 1 




/X+y 


(X-y) 


= 0-i)Tr'/^ n-O+i) E* [2y(2-x)]} 


(2.29) ; 


where E* is the complex conjugate of E as defined in (2.21). 

(2) 

Substituting for l(y,x,0 + ) in p* ‘\x<2,0 + ), we have 
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I 


I 


j 


t 


I 


I 


yiV"^ 

from which (2.25) becomes 


{l-(VH)E*r.2ii(2-x)]> (2.30); 



so that the Kutta condition is satisfied at x=2 as required. 
The result in (2. 3D agrees with that given in ref. 13 for 
the pressure distribution on an infinite-span wing passing 
through an oblique short-wavelength gust. It serves as a 
check for the expression in (2.24) that we have derived for 
the pressure p(x-,Z) at a general field point. 


Two-dimensional acoustic far field in the airfoil frame of 
reference 

Here we start with the two-dimensional pressure field 
P(x,y,Z,t) due to the passage of an infinite-span wing 
through an oblique gust as given by (2.24) and (2.13a). We 
determine the two-dimensional far field in the airfoil 
reference frame by allowing the field point to go to 
infinity. The calculated directivity of , sound is 

theoretically that which would be measured in a wind tunnel 
experiment. 

From (2.24) and (2.13a) we have 


I 
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p (x-,y.Z,t) = - PnWp u z /lzl 
2 tt /TW /'i^ 


iojt - ikyy + ikxM^x 
• ^ V-M^ 


V 





1 e 


2iX 


C ^/X+u •’y V c - y 

We notice that in order to obtain a nonvanishing far 


Jy ^ ^ 1 . 




field the largest contribution to P(x,y,Z,t) as/x^+Z^-m* must 


come from values of X such that/X*^ <o or |xl '<y, and so it 
follows that in the far field, the largest contribution to 
the inside integral in (2.32) comes from its lower limit. 


Integrating by parts, we thus obtain that as/x^+Z^^oo 


C-X 




]i 


c.-y 


= TTv^ e'^''^{l-(l+-i).E*C2(y-X)]} (2.33) 
y-x 


Substituting this into (2.32.), the asymptotic behavior 
of the X integral as /x^ may be determined by steepest 


descents (Noble, 


(23) 


pp. 33-36) to be 


2 ^ ^/iW- jlilTizyF- 




1 - cos (.0*/2) - (1 + i) E* [2y (l-cose*)]") 


0 - M2 sin20 )'/■* 




where the substitution Z= A -M has been made so that 

_ 1 I ( ■ " 

9 = tan (z/x) and 0* =tan ( /I -M^tan 0) . The expression in (2.34) 
shows that for negative y^ (y =-i/- y"") there is no acoustic 
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field (the solution in (2.34) is exponentially damped). We 
had anticipated that this should be the case since the 
disturbance due to the passage of the infinite-span wing 
through the gust then travels subsonically through the fluid. 
The quantity in brackets is the two-dimensional acoustic 
directivity D(Q) in the airfoil reference frame; we plot it 
for different values of the Mach number M in figs. 2.5a, b. 
For values of 1 foot for the semichord b and 1100 ft/sec for 
the sound speed Cq, the acoustic signals corresponding to 
these patterns are sound tones of 280 and 315 Hz, 
respectively. 

We can appreciate the dramatic effect of noncompactness 
on acoustic radiation if while we look at these we recall the 
directivity pattern for a simple dipole - the acoustic field 
for a compact region. Note the sensitivity of the lobing to 
a relatively small change in Mach number. The parameter 
which actually controls the amount of noncompactness of 
course is u > which lumps frequency and ‘compressibility 
effects together and without which the Helmholtz equation 
(2.3) becomes the Laplacian, the governing equation for 
compact aerodynamic regions; it appears alone in the 
directivity part of (2.34). 

(j 6 ) ( I 3 ) 

Amiet and Adamczyk' have compared the predicted 
aerodynamic loading (2.31) to numerical results obtained by 

(2i») 

Graham for a wide range of y values. They have determined 
that the two-term theory (one leading- and trailing-edge 
problem) gives good results for loading when y >.7, 



approximately, for y real. For imaginary U (2.31) remains 
valid in a wider range, apparently because the trailing-edge 
part of the loading converges then much faster; the Fresnel 
integral turns into an Error function and the loading- becomes 
hydrodynamic in character. 
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Chapter 3:. Approximate aeroacoustic model for a rotating 
blade of finite span passing through a gust. Application for 
a potential vortex 

In chapter 2 we derived an expression for the pressure 
field everywhere in space for the two-dimensional interaction 
of a subsonic infinite-span airfoil with an oblique gust. In 
this chapter we use that two-dimensional solution in a 
superposition procedure which will enable us to model the 
three-dimensional acoustic field radiated by a rectangular 
blade rotating with high subsonic tip-Mach number. We find 
that when the disturbance due to the interaction travels 
supersonically through the fluid, the model breaks down along 
certain directions. For cases of subsonic trace speed of the 
disturbance, we present typical examples of acoustic 
directivity patterns and compare one case to the 

corresponding two-dimensional result in chapter 2. 

(IS ) 

Amiet has pointed out that the contribution to the 
acoustic field from the acceleration effect of sources of 
sound on a rotating blade may be neglected when o) , the 
acoustic frequency, is much larger than the blade angular 
speed. Our three-dimensional model for the acoustic field 
uses this simplifying fact, and so we assume that each 
spanwise section of the blade passes rectilinearly through 
the gust with local speed varying from hub to tip. 

In addition, the associated three-dimensional 
aerodynamic model makes the following two related 
assumptions; (1) the chordwise unsteadiness dominates the 
three-dimensional tip effect over all of the span; and (2) 
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that the spanwihe loading may adequately modelled by a 
triangular shape with maximum loading at the tip (in reality 
the actual loading near the tip must vanish as y'L-y ) . 

We close the chapter by applying the leading-edge part 
of our result to predict acoustic signatures for the 
interaction with a potential vortex. The vortex in the model 
is located a distance h chords below the flight plane. 

We start by defining the three-dimensional pressure 
field jj(x,y,Z,t) clue to the passage of a rotating blade of 
span L through an oblique gust; as before Z= /I z. We 
require that this three-dimensional pressure equal the 
two-dimensional pressure P given in (2.32) on Z=0+, y=L, and 
that it decrease from this value to zero linearly at Z=0+, 
y=0 . We make the changes 


^3-D(x,y,Z»t) = P* (x,y,Z) 6^“"^ 


and Ysvf-M^y so that P* satisfies 


^*xx ^*YY = 0 


subject to the boundary condition 


(3.1), 


(3.2) 


P*(x,Y',0+) =( J_ 


'L/MF 
0 


P* 2 _o (x.Y,0+) for 0<Y</nF L 
otherwise 


(3.3), 


where 


P* 2 — Q 0+) — ~PnWo U 1 s 


Ax >... 


• ikx tan A Y 
/wF 
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A-y uA!+y u ^ 
as may be seen from eqn. (2.32). 

We now define the transform pair 


dC /l^ e 

C-x V c-y 


-215 


p*(x.y,z) = 


f dX,dX2 fiXuXiil) 

(2u)2 


P- Ui.XzjZ) = 


cjx dY e 


iXiX + iXaY 


P*(x,Y,Z) 


Transforming (3.2) and solving, we obtain that 


P CAi,A 2;Z+)= + P (Xi,X2,0+) ZAf +Xii-K" 


where we have called kj^M/( 1-M ^IbK, for convenience. 
(3.5b) and (3.4) it follows that 

P (Ai,X 2,0+) = F (Xi) Q 1 X 2 - tanA ^ 

1. /wF"j 


where 




/TFFl 


A 2 - ky tanA 

/MF 






V 


and 

F Ui) 


~ - [J_ 


/PF /t: 


1 e 


2iXi 


ia^y (j/Xi-y irAi+y 


/€+ y e 

C-Xi V y 


y 


(3.4) . 


(3.5a), 


(3.5b). 


(3.6) , 
From 


(3.7) , 


(3.8), 



(3.9). 


From (3.7), (3.6), (3.4a), and (3.1) we have 


theref ore 
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that 


„ / V -7 j \ y>it ^ ikyH^x 

P3_D(^»Y>Z,t) “ G _ “T^FP” 

dAidX,. 


(27t) = 


F (Xi) G(X,- kytanA j 

/TW I (3.10) . 


We now translate the solution in (3.10) from the airfoil 
frame of reference to the ground, or acoustic, reference 
frame. We let Xi=bXf-l<^ / ( 1 /1-M^ , and express the 

result in (3.10) in terms of the original unstretched y, z 
variables. We also introduce a oj' -t Fourier transform so 
that 

, fOO I « > 

^ .giVt-iX.bx-iXaby 




L 


^d5:id^2 

(2.-ny 


exp^ibz 


OJ 


12 


I 


-X1-X2 


5 [aj'-(a>-^iU)] F (b3(i+ki NP) G 

T-FP“ 


bXg-ky tanA 


/TFF 


(3.11) , 

where k'x=(oj+XjU)b/U. Next, we express the 6 -function 
appearing in (3.11) in terms of one of its well known 
limiting forms, e.g., 5 (a)=lim 1 -jrrr and interchange the 

orders of XyXa integration with the u)' integral: 

CO 

" s I j. r r . . L 

P 3 _o(>J^>y>- 2 .t) = Lim ^ 


£-^■0 (2Ttf 


-ibz 


r^ji2 ^ 


• , - 103 ' t f 
QUJ G ' 




1? TF?nr^ls:ix;i)T>} 


bXz- ky tanA 


/i -M 2 

(3.12). 

Here bx_=bx-Ut, the "a" subscript denoting the acoustic 

a 


reference frame. The last step is to evaluate P 


3-0 


in the far 
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field; this may be readily done by approximating the Xj, \z 

(2S) 

integrals using the method of stationary phase (V/hitham, p. 

392). Before we do this, however, we may use the 

approximation indicated in (2.33) for the integral appearing 
in F in (3.9). The final result for the sound signal in the 
three-dimensional acoustic field which a listener on the 
ground hears as the rotating blade of span L and tip Mach 
number M passes through a stationary gust of short wavelength 
is 



(3.14a) , 


and 



$ 
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y/r^, z/rg specify the position of the observer on the 
ground. D (x^/r^, y/r^, z/rg)' is the three-dimensional 
acoustic directivity not including the factor (1+Mxa/ra) 
which gives the usual forward enhancement of the acoustic 
signal from a moving source, in this case, the blade 
traveling in the -x^ direction; l+Mx^/rg also represents the 
position-dependent Doppler shift of the tone, as indicated in 
the exponential in (3.13). 

The expressions in (3.13)j (3.14a) stand as are for 

M/sinA>1 ( tanA<M/ A ) , that is, for supersonic trace speeds 
of the disturbance. For subsonic trace speeds, sin >M, and 
we recall from chapter 2 that u =-i/-u^ , so that may then 
write 



in (3.13), (3.14a). Before we make this change for the 
subsonic cases, however, we notice that the fiVst term inside 
brackets in (3.14a), the contribution from the leading-edge 
problem, blows up in the direction given by 
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As fig. 3.1 shows, this is the direction taken by 
acoustic rays born at the line formed by the intersection of 
the downstream part of a typical Mach cone and the x^-y 
plane. The breakdown of the leading-edge solution in this 
direction can be explained mathematically as follows: Since 
the chord extends downstream to infinity, in the' above 
operations we have in fact tried to calculate the acoustic 
field of a. distribution of supersonic sources which is 
infinite in extent; the surfaces of the associated system of 
Mach cones do not interfere and so, even those farthest from 
the leading edge will have an acoustic effect. But the 
determination of such an acoustic field is a badly posed 
problem for the method of stationary phase, which relies on 
the existence of a single point in physical space (and also 
in wavenumber space) where the acoustic effect of the entire 
source distribution may be modelled by one equivalent source 
of arbitrary complexity. 

Similarly, we notice that the second term in brackets 
in (3.14a), the contribution from the trailing-edge problem, 
blows up along two different directions: one of these is 
given in (3.15), and the other by the solution of 



ra 


which is 
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They are the directions of acoustic propagation from the two 
lines formed by the intersection of the Mach cones and the 
Xg-y plane. That the trailing-edge part of the solution 

breaks down in both directions should not surprise us if we 
recall that although it was constructed to cancel the 

leading-edge result in the wake, it was also nonzero (though 

small) upstream of the plate. Since the trailing-edge 

I solution cancels the leading-edge solution in the wake, 

however, we would expect it to cancel also the breakdown 
because, in effect, it changes the leading edge’s region of 
1 semi-infinite extent - the blade chord plus the wake, into 

one of finite extent - just the blade chord. Below we verify 

that this is in fact the case. Unfortunately, we still end 
’ up with a part of the trailing-edge solution which breaks 

down in the upstream direction given by (3.16); to correct 
it, a second leading-edge problem would have to be solved. 
But even then our problems would not be over, for we would 
still have a semi-infinite region of supersonic sources, now 
in the downstream direction. 


In the direction given by (3.15), we find that 




4.3 


so that the part in curly brackets of the trailing-edge 
contribution to the directivity pattern in (3.14a) (the 
Fresnel integral vanishes since its argument vanishes) 
becomes 


VPFP-. " 


(1 + M ^a) 


1-M^ 


2 

tanA 


1/2 


thus cancelling that from the leading edge. If we now look 
at the high-frequency behavior of the term inside the large 
curly brackets, we obtain that for large k , it becomes 

A 

1 


rri 


/2 I MV . 2 .‘ ^ 
- tan 2 A 


, n ^ 


\ 


® 




- tan^A 






(^4 * M) 

/r^ 0 + M xa/ra) 


1/2 


VJe conclude that in the high-frequency limit k^^soo, the 
far-field we have calculated for supersonic trace speeds is 
well behaved everywhere, even along the direction given by 
(3.16). However, since the limiting processes cannot be 
interchanged, that is, the solution is not well behaved if we 
first approach the direction of (3.16) and then let k^ , we 
note that for supersonic trace speeds the uniform convergence 
promised by Landahl^^^^ in two dimensions does not carry to 


three. 


We notice from (3.13) and (3.14a,b), that/for large k^ 
the predicted acoustic spectrum (not including w^ ) decays as 
1/k;^., the behavior of the second term on the right side of 
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(3.14b) in the expression for . this large-k^ behavior 
has an associated small-time behavior for the transient 
signal which may or may not be correct - recall that the 
present aerodynamic model does not predict well radiation 
from tip dipoles. However, as we will see in chapter 5, 
where we develop a better aerodynamic model for the blade tip 
edge, the l/k^ behavior is the correct one after all; so that 
in fact, for small times the expressions in (3.13), (3.14a,b) 
are off by at most a multiplicative constant. 

Since the trailing-edge part of the solution hardly 
contributes for small times (the high-frequency beh-avior is 
given by the leading-edge solution, whose directivity does 
not depend on frequency), essentially all the phase 
information in contained in GZ/I-M^, which can be traced back 
to the I 2 transform of the spanwise linear loading. Also, 
since G/ /! -M^ =i/2 at k =0, we see that the predicted 
spectrum takes on a constant value at 

We now express (3.14a, b) for the case of a subsonic 
trace speed of disturbances. Also, we introduce the 
spherical coordinate system (r, ,y, 4 )) shown in fig. 3.2, in 
order to define the position of a listener on the ground. 


We have then, t^at 




- € - 




E* 


/T-M^ 



tan^A- 


_M (sincfccosrH'l) 

U+M sin(j)Cosv) 


Msin<t)sinv 
1+M sin4>cosY 


- tanA 


//I 



/Msin^siny 

Ky 1 I m i J -'T'i 


l-H-1 sin4)COSY 


tanAl 


//TmT - -iXL (tanA- Msind>siny ) 
//i-M - _ X l+flsin(J)COSY 


-1 


. -ikxL (tanA- Hsinc&sinY 
+ ^ I+Msintbcosy 

kw (tan A- Msin4>sinY/(l-H-lsin(i!COSy) 


kx [tanA - MsinctsinY 1 
t l+Msin(j)COSY'^ 


(3.17b). 


For large A, the imaginary term in the argument of the 
Fresnel integral dominates the real part. As a result, the 
Fresnel integral has essentially an Error-function behavior 
which, as a function of y , does not provide for much lobing 
in the directivity pattern - certainly not as much as we had 
in the two-dimensional patterns of chapter 2. This .is 
generally apparent in part (a) of figs. 3. 3-3. 8 for the 
three-dimensional directivity function D 3 _jj(Y , li*) for a number 
of cases of interest. Actually, in order to make comparisons 
with the results from chapter 2, we have plotted here the 
quantity 

^ D3-0 (Y»j>) 

2 /ir (l+Msin({)COSY)^ 

The difference in the acoustic spreading effect between two 
and three dimensions (the ^/ /F and 1/r, factors) is not taken 

a 

into account in the comparison. In order to include this 



effect, we would multiply the above expression by 

In fig. 3.3a-d we have selected a set of physical 
parameters (M,A,oib/U) so that the resulting value of Ip^I 
matches that of fig. 2.5a. Notice the decreased efficiency 
of the wing as an acoustic radiator for the present 
three-dimensional case with load fronts which travel 
subsonically through the still fluid. The cos ^ G/»^ 1-M^ part 
is then a slowly varying function ofY,4>; essentially it has 
a radiation pattern on the plane Y =90® , -TT/2<<f><n/2 (normal to 
the flight direction) equivalent to that for a strong tip 
dipole plus a somewhat weaker dipole located at some distance 
inboard. On the plane parallel to the flight direction 
(y=0®, -tt/ 2<4) <ir/2) , the contribution of cos ({> G//T^M^ is just 
that of a simple dipole. As previously discussed, for the 
two-dimensional case from chapter 2 the speed of the 
disturbance through the fluid was supersonic, and so, every 
section along the blade's infinite span radiated noise. 

Figs. 3.3, 3.4, and 3.5 show the differences in the 
directivity patterns for three different tones. The 
"rippling" increases with reduced frequency, and the 
magnitudes decrease. Figs. 3.4 and 3.6 illustrate the effect 
of a change in Mach number. Compressibility has the same 
effect of increasing the rippling of the pattern - again, 
bringing out the lumped character of compressibility and 
frequency effects for our noncompact blade. Another 
important high-speed effect is the increased forward 
enhancement ofsignals (compare fig. 3.4a to fig. 3.6a), 
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From figs. 3.5 and 3.7, we see the effect of a change in 
A, the interaction angle, on the radiation field. Notice 
that the greatest qualitative differences occur on the 
azimuth y=90°. This is consistent with the . earlier 
interpretation of the radiation field on the plane normal to 
the flight direction as essentially due to a strong tip 
dipole and a weak inboard dipole: changes in mostly affect 
their difference in phase. Also, notice the reduced levels 
for a=809 The reason for this is somewhat more subtle: for 
larger , the effective trace Mach number of the load fronts 
has a lower value; as a result G/ lA in (3.17b) is farther 
for A=80° than for A=62° from assuming its peaky maximum 
value for a given frequency. Physically, there can be more 
cancellation of acoustic signals for the greater value of 
A=80^ 

Finally, in fig. 3.8 we consider a case in which the 
trace speed of sources through the fluid is nearly transonic: 
we choose our parameters M,A,ojb/U so that takes on the 
relatively small value of -1 . Since the imaginary part of 
the Fresnel integral in (3.17a) is then small ( tanA=M//T=TT*", 
or sinA=M), we would expect the latter to behave as one, 
rather than as an Error function. The observed complexity of 
the pattern for y=l80°, 0<<P<tt/2 clearly shows the increased 
participation of the trailing-edge part of the solution. For 
the y=135°, 150° azimuths, we notice that G/ is reaching 
its peaky behavior, indicating physically that a preferred 

direction of strong radiation is about to appear normal to 
the load fronts (the transonic Mach plane). 
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Application of leading-edge theory to prediction of noise for 
a potential vortex 

We consider a potential vortex in the still fluid 
oriented at interaction angle A, measured from the leading 
edge (fig. 3.9a). The vortex is located a distance h chords 
below the plane of flight, where it induces the upwash 
pattern in a plane normal to its own axis indicated in fig. 
3.9b, and given by the Biot-Savart law; 

w(C)=-l7b g 
2rr 

(3.18), 

where V is the vortex strength, in ft^/sec. By Fourier's 
theorem, we may express this upwash as a sum of sinusoidal 
gusts 


w 


(g) = 


dke 


-ikc 


W(k) 


(3.19) . 


Since ^=XgCosA+ysin A, we have that 
W(k) = 


J 

00 


f to* 

(2u] 2 v/(Xj^C0SA+ysinA) =S (ky-kxtanA) 1^(kx/cosA) 

cosA 

( 3 . 20 ) , 


from which 


w(xaCosA+ysinA) = 


dkydkx 


e"ikx-Xa"^*^yy 
W (k) = 


® W Ux/cosa) 

co^A (3.21). 


If we now recall that bx^=bx-Ut, where x -is' fixed on the 
wing, we see that (3.21) becomes 


I 
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w(xacosA+ysinA) 


00 



dkxe 


ikyUt 


- ik^x - iK tanA y 
Vq Ocx/cosaO 


(3.22), 


where we have called W(kx/cosA)/cosA =w^ (kx/cosA ) , and 
identified it as the magnitude of the convected gust .(as seen 
by an observer on the airfoil) in all our work up to this 
point. Notice that k^U/b=o), the acoustic frequency if 
kx=ccb/U, the reduced frequency. 

In chapter 2 and the first part of the present chapter 
we have determined the tone emitted by the passage of a blade 
through a gust, the integrand in (3.22). In order to 
calculate the acoustic field for the passage over the vortex,, 
we simply superpose the solution of the gust problems 
properly weighted by the function WQ(kx/cosA). We have then, 
for the transient signal ptr,t), that 


P (r,t) = 


CO 

dkx P(r,t;kx) 

I 


(3.23). 


We now want to evaluate p(r, t) for the model developed 
in this chapter, that is, for P(F, t;k^) given by (3.16), 
(3, 17a,b). We have shown that for the high part of the 
spectrum P may be adequately approximated by the leading-edge 
part of the solution, which we rewrite here for convenience 
with a few changes in notation: 


P(r,t;k }= - ipn U 
2tt 


Wo(kv/cosA) /A-cosa 

y^l-i|tanA| >1 ^p(l+MsinqcosY)^ 

^'sTtFa j 
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Ihxl. 

e b 


^ lltanA^ 

I 

(' 


172 


w 

(l+MsJntjjcosY) 


dy' e-% 

r 


tan A- Msin4>cosY \ 
HMsTn4>cosY/ 


(3.24) , 


where r=t-b^/cQ, the source, or retarted time. 

If we now let x* =t/ ( 1+Msin cfcos y) » substitute for P into 
(3.23) and apply the convolution theorem, we find that 


p(r,T*l = 


,^g(T') f(-T*-r') 


2'tt 


where 


g(.T*). = 2-rr* for o<x*<L tanA- Msindsiny 
L l+Msin(j>cosY 

CO 


fCx* I = 


Jij _ ikxT* 

% Wq (k^/cosA) 


For a potential vortex, we find that 


% Cgiil = - wu ^ 


4ttcosA kv 


zh |k 

cosA ^ 


so that 


(3.25), 


(3.26a) , 


(3.26b). 


(3.27); 


f(x*)= r/b X* ^ 

ZttcosA t* 2+ (h/cosA)^ 

K 5,do ) , 

With f and g so given, the integral in (3.25) may be 
evaluated in closed form. The final result for the acoustic 
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signature which a listener on the ground hears bi|/Cg time 
units after the blade has passed over the vortex is 


p (r,T)=-i Po ^ ( 1 ^ 

2u 


1 -i 

itan^A- ) 


1 w 3 


,M(sin<t>cosrtMl + iitanAl /TW \ ' 

»UMsin^cosT 


r/b 

/ ^ LcosA 


T_ 

ra 




COS(t> 


(l+Msin4)COSY)* 


tanA- Msinc^siny \ 

l+Msin(|)COSY ] 


..h arctan 

Ty{l+Msin<l)COSY) 

h arctan 

r T 1 

{tanA-Msi nd)Si n y \ 

cosA 

h/cosA 

cos A 

l+Msin^jCOSY 

h/cosA 

‘ l+MsinicosA/ 


-■^/C 2 (l+Msin<J>cosY)J 


log 


l+f^in<j)COSY 


- L 


ftanA-M sindiSinv) 

1+Msin(j)C0SY 


(h/cosA) 


(»-M ) 


31. r 


(3.29). 


£f/(l+Msin({)COSY)]^ + (h/cosA 

X — "CLK/ia 

For the present noncorapact situation, there is no 
guarantee that the maximum values of the acoustic signal 
(3.29) will coincide with the field point where its tones 
individually reach their maximum amplitudes. From figs, 

3. 3-3. 5, we see that the far-field point where each tone had 
its maximum is not a very strong function of frequency; so 
that essentially they all had large amplitudes at 
(f=65® , Y=1 80®-, for the cases in figs. 3. 3-3. 5. Still, the 
phases of the spectral components could conspire to cause a 
great deal of cancellation, and the far-field point of 
maximum peak level might be where we least expect it; at some 
point where the tone amplitudes are not particularly high. 
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but where their phases reinforce strongly. 

Such fears and suspicions prompted the calculation of 
part the directivity pattern, shown in fig. 3.10, of the peak 
signal for the case of M=.8, A=62°, the case corresponding to 
the tones of figs. 3. 3-3*5. We used a value of 238 ft Vsec 
for the maximum sectional circulation (assumed occurring at 
the 75? span). For M=.8 this choice implies a maximum 
sectional lift L. p^c^r = 498 Ib/ft, which, if assumed as 
acting at the blade tip and decreasing linearly to zero at 
the hub, gives a value of 3735 lb of lift for each 15-ft 
blade. 

We find that although the point <jp6? ,y=l80® is in fact 
where the signature reaches its maximum level, substantial 
variations exist: for example, the peak level at (j)=70® ,y = 1 5(f 
was only 114 dB for h/b=,1, compared to that at ({>=65° , y=1 80° , 
where the signal had the maximum level of 131 dB for the same 
value of blade-vortex separation h. From just looking at the 
patterns in parts (a) and (e), for (p=65°, y=70®, 
respectively, of figs, 3*3-3. 5, one should never have 
expected such a large difference. 

Fig, 3*11a,b shows predicted signatures at the point 
(j)=65® ,y = 1 80° , where the signal was loudest for h/b=.1 (fig. 
3.11b). The value of r is 1000 ft. The case for h/b=.5 
(fig. 3.11a) had a lower peak value of 126 dB. This is to be 
expected, since dF/dt on the blade should then be less. 

The purpose of fig. 3.11c is to show, upon comparison to 
fig. 3.11b, that our closed-form solution (3.29) contains 


1 
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the Doppler-shif t information of the blade-slap signature. 
As indicated, the signal in fig. 3.11c was evaluated at a 
point also on the plane parallel to the flight direction (for 
Y=180*^, 0 <'K‘it/ 2), 'but a~t a higher elevation angle so that 
it should be less shifted. And it is. 

Fig. 3.12 shows the spectra of signals corresponding to 
the M, A cases of figs. 3. 3-3. 8. They were calculated at the 
field position where each tone reached its maximum amplitude. 
As previously noted, for M=0.8, ^= 62 ° this point coincided 
with the maximum peak level; the same may not be true of the 
others. The main objective of the figure is to show the 
large differences which exist in the rates of decay of the 
spectra between the h/b=.5 and h/b=.1 cases, and to indicate, 
for one of the cases, the value of frequency above which the 


present noncompact model 

is 

valid. 

We 

recall from 

the 

discussion at the end 

of 

chapter 

2 

that effects 

of 


noncompactness become significant for y >.7, for supersonic 
trace speeds, and for u- greater than some value less than 
.7 for subsonic trace speeds. For A=62® M=.8, P*>.7 for 
f>8l9 Hz. This is the section of the spectrum which should 
be predicted using a noncorapact theory; as we clearly see, 
for h/b=,1 it contains a substantial amount of the 
signature’s acoustic energy. 



Chapter 4; Quarter-infinite plate with side edge passing 
through a gust 

The acoustic model developed in chapter 3 overestimated 
the strength and phase of acoustic dipoles near the tip of a 
rotating blade. This must be so, since there we simply cut 
off the loading abruptly at the tip, where in reality we know 
it must vanish as vL-y. Yet, when the load fronts due to the 
interaction travel subsonically through the fluid, that is, 
for M/sinA<1, only outboard sections actually radiate sound. 
Because this type of interaction is generally believed to 
occur often in single-rotor helicopters, it is important, 
therefore, that we be able to model the strength of tip 
dipoles more accurately than we did in chapter 3. That is 
the purpose of the present chapter. 

The model we develop here is more rudimentary in some 
ways than that worked out earlier: for simplicity we remove 
the trailing edge and allow the wing to have a semi-infinite 
span; also we cannot model blade rotation here. Still, the 
theory we develop below for a quater-inf inite plate passing 
through a gust proves itself useful in chapter 5, where we 
make a rational simplification which allows us to reintroduce 
the trailing edge and so, to predict acoustic radiation from 
a blade of semi-infinite span with correct side-edge loading 
passing through a gust. 

We start here by stating the boundary-value problem for 
the three-dimensional velocity potential fiel-d-'-$, due to the 
interaction of a quarter-infinite plate with side edge with 
an oblique gust. Again, the oblique gust is that wavenumber 



55 


component of amplitude Wq of the upwash on the plane of 
flight induced by a vortex convected at angle of obliqueness 
A by a subsonic freestream (fig. 4.1). The governing 
equation - the linearized convected-wave equation, and 
the boundary conditions are 


(J) +$ + A 

""xx ^yy ^Z2 






-WoS y 


(4.1), 
(4.2a) , 


^Cx,y<0,0,t) = $(x<0,y,0,t) = 0 


(4.2b); 


where we use the symbol for the reduced frequency as we 
did in chapters 2 and 3. Here however, k^=ujs/U, s being a 
characteristic length in the problem which at this point we 
leave unspecified. The reason for this ambiguity is that 
although the only characteristic length for now is the 
acoustic wavelength itself, in chapter 5 we ’ reintroduce the 
trailing edge, and with it, the blade semichord - a more 
useful characteristic length for the model. 

The three-dimensional pressure field P(x,y,z,t) is also 
required to satisfy (4.1) and (4.2b). In addition, we expect 
it to have the following behavior near the edges 

P((-^0+,y>0,0,t) 'V 1//X (4.3a), 


P(x •>0,y->-0+,0 ,t)'vy 


(4.3b) . 
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As before in chapter 3, we let $ = (^xpCiCcu t+k^^M ^x/( 1-M^')] ; 
notice, however, that the exp[-ik tanAy] term is not 

X 

included. V/e make the changes y=Y/ /T=7F, z=Z/i^1-M^ . 

Similar changes are made for P. The boundary-value problem 
fortj)* becomes 


*^*xx ^ ° (4.4) , 

with 

<(>*_(X>0,Y>0,0) =-wo$ (4.5a), 

<{)* (x<0,y,0) = <j>* fx,y<0,0)= 0 (4.5b). 


I 


We define the transform pair 

dAi e <^(Ai,A 2 'j 2 ) 

(4.6a), 
Cl 

<P,(xJ,Z) (II. 6b); 

where Ci, and Cc denote complex integration contours in the 
Xi , X 2 planes, respectively. We apply (4.6b) to (4.4). The 
solution for positive Z is 


4>*(x,Y,Z) = 




dx, e 


-iX2Y 


4>(A,A2’>Z) = 


dx dY e 
Ztt 


iA 


d Ai e 

/2T 


iAi X 


d A 2 
/ 2tt 


•iA^ - ZAt+A!-K" ~ ^2;0+)(4.7) 


4>*(x,Y.Z) = 
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Since (jjit(x,Y<0,0)=0, it follows from (4.7) with Z=0 that 
<p(\nX^ must be an analytic function of X 2 in the upper half x '2 
plane. Using standard nomenclature, we therefore say that 
4 )=(j) 0 . From (4.7) we may calculate |^ ( x , Y, 0 ) , the vertical 
velocity field from the solution on the plane of the plate; 

(x,Y,0+)=-f ^ dX 2 e"''^ 2 Y ^Xi,X2;0+) (4.8) . 

Z /-s- ,-5- (±) 

r» / 2u -/ 2tt 

I 

The boundary-value problem posed by (4,4), (4.5a,b) 

cannot be solved exactly for subsonic flow by any simple 
analytical procedure known. The assumptions which follow, 
however, allow us to obtain an approximate solution which 
later we show is valid and useful for short acoustic 
wavelengths; effectively, by making them we relax the 
requirement of (x<0, Y, Z) =0. Vie let 

4 (Xi,X2,0+) AT+Xpl^ = r(Xi,X2) + v(Xi,X2) (4.9), 


where 



r(x,Y) = 
v(x,Y) = 


0 for Y<0 

2-0 UPWASH for Y>0 


? for Y<0 
o' for Y>0 


C4.10a) , 
(4.10b). 


From chapter 2, equation (2.12), the upwash on the plane 
Z=0 of the two-dimensional solution turns out to be 
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r(x,Y,0) = -iwns e 


-ikxtanAY 




/5 

V 1- 


,M2 + w 


dXi 

(Ai- kx 0 


-iXx 

® . >/x7nr(4.ii) , 


Ci 


where Cj is as C in fig. 2.2. We therefore find that 


r(Xj ,X2) “ A(Xi) 


Xz- kytanA 


(4.12), 


where 


A(Xi) = - WnS/Xi-i-y 


2-ir/pF 


1 

(Xi- kx ) 


In the solution we are constructing, v(x,Y>0)=0 as 
stated in (4.10b). It follows, therefore, that v(Xi,X^ must 
be an analytic function of X 2 in the lower half X 2 plane, and 
so we write v(Xi,X 2 ) =Vq(Xi,X 2 ) . Substituting for r and v in 


(4.9), and factoring the term /Xt+Xt-K'^ , we have that 


<j>(X2)A“2-(K2-X?)V2= A(Xi) 
© - • 


A(Xi) 1 I UX.) 

X2-Mm\ A 2 + [K'^-X'iPJ^ /X^'(K^i-xi)'A ' 

where we have written d>^(X^ instead of (b^XuXz^ to avoid 

© © 

confusion as to the meaning of the symbol 0, which refers 
here to the analyticity of (p in the upper half of the X 2 
plane, not the Xi plane. 

We recall now that since an exp(iwt) behavior for the 
harmonic time dependence was made, the condition of outward 
wave propagation requires us, to allow m, and so ky and y, to 
have a small negative imaginary part. Also, we notice from 
( 4 . 13 ) that the two radicals /X 2 -(K'-Xtj ^7 and /X 2 +(K‘^-Xi )^'"contain 


I 
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the quantity , a function of Ai which has branch points 
at Aiz±K. We choose that branch of ATTfp’ with vanishing 
argument as along the Ci contour. This choice now 
implies two new branch cuts in the \i plane, which now is as 
shown in fig. 4.2a. Fig. 4.2b indicates the position of the 
singularities in the X 2 plane and the corresponding regions of 
0 and 0 analyticity. 

If we split in the usual way the first term on the right 
side of (4.13) into a sum of two functions, one 0 and the 
other 0 , we obtain a Wiener-Hopf equation of the form 0 = 0 
from which we may find a straightf orward manner; the 
details parallel those for the leading-edge problem in 
chapter 2 and so will not be repeated here. We find that 



Substituting for A and by means of (4.7) we obtain 
(^i(x,Y,Z). From the linearized relation between pressure and 
velocity potential p*=-(PqU/s ) ( ( ik ^/d )<^^+9<fiy3 x) , we may 
also calculate p^(x;Y,Z+). The final results are; 
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!>-*■ . (x,Y,Z) = - ipnWn u_ 1 

/jr 




+ u 


-Uix ^ 

dXie yXi+ij 


Cl 


i^xtanA/ + ( i^2_;^2 j 


TJT 


(A 2 - k^tanA/Zi-M^.) 


V 


(4.16). 


/xT^F^irr/^' 

Asymptotic behavior of solution near the edges 

By construction, the integrand of the X 2 integral is an 
analytic function of X 2 in the upper half X 2 plane, and 
consequently p*(x, Y<0,0)=0. However, the integrand of the Xi 
integral is not an analytic function of Xi in the upper half 
Xi plane. As a result, the solution in (4.16) does not 
satisfy the condition of p*(x<0,Y,0)=0. The next step in our 
analysis is to determine how severely the result in (4.16) 
violates this upstream condition and also to find whether it 
has the required behavior near the edges. V/e therefore set 
Z=0+ and obtain from (4.16), after evaluating the X 2 integral 


-ik^._tanAV 


in closed form, that 
P *(x ,1V 0,0+) = (i-1) rhWn U e 


-Uix. 


■ixx , 


+ u 


dXi R K 


v'X 1 - y 


Y 

L V 


fkji^an^ 




-M" 


/R"“xr 

(4.17) 


In order to investigate the solution for small values of 
X, that is, near the leading edge, we look at the large-Xi 
asymptotic form for the integrand (at the same time with Y 
not small). Expanding the Fresnel integral for large 
argument, we have 


p*(x,Y>0,0) = - PoWn y_ 
2ir 


1 




/! 




-ikvtan Y 




+U 




.. dXi e 
/Xr U 


-iXiX 


t 
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■e . V (-i) (2n+l)! 1 

;^0 W»T7i-| /ahT ^ 

so that the leading term in the behavior near the leading 
edge is the two-dimensional leading-edge solution in chapter 
2's equation (2.15), which vanishes for x<0 and behaves as 
1/*^ for x=0+. The part of the solution which does not 
vanish upstream of the plate appears as an asymptotic series 
in inverse powers of Y, the distance from the side edge 
normalized by the characteristic length s; it also represents 
the three dimensional effect of the side edge on the loading 
p*(x,Y>0,0). 

To calculate the coefficients of the asymptotic series 
in (4,18), we must evaluate the integral 


In(x,Y>0) = 


-iXiX-Y/xfF 
dX, e _J 

(/TW J 


(4.19) 


for Y>0, x->0+. By the convolution theorem, we may rewrite 


it as 


00 

I„(^.RO) = - 1 /? 

, 


(K ^ 


-fn^5)(4.20) , 


where 


dXj 1 

^ (4.21). 

//HF • V 


The integrand in f^ (x) has branch-point singularities at 


X 
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Ai =+K. Besides these, we also find that as Xi-^±y 


1 


kxtanA , 


^n+i /2 ^ 2kytar.A _,w, 1 , 


so that in addition we have a pole of order n+1 at Xi = y and 
another branch point atXi=-y. We write therefore that 


= f^^x) + ff (X) 
!„(*) = li”(x) + li'^(x) 


(4.23a), 


(4.23b) ; 


where 


fi“(x) = ^ 




(X) - ^ 


^ 2 kytanA 
/l 




and 


Branch line 
at Ai=+K 


U\x fkytanA T^TT^n-) 

1 1 1 np V J 


(4.24b) ; 


I^(x,V) = -1£ 


% 


.( 2 ) 


d^Hr^ IK/ Cx-Cj^+V^) . KY 


fi^(? ) (4.25a), 


If'^(x.y) =-i y| 


— CO 

fOO 


d^nS^^ (K/ (x-c)^+ Y^) KY 


f^'^(^)(4.25b). 


- oo 

■Jy 


/U^+ 

For n=0,1 f|^(x) can be calculated by deforming the 
contour Cj in fig. 4.3 to C* , or , for x>0, x<0, 

i ■ 

respectively. 


We then obtain 
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1 a. . I / o M -lUX \, iur A-'u 

<.+u 1 /■ZkytanA ntt {-2:n) e 4 ! 

W ’ 'JE b - E*[2 u(-x)] for x<0 






for x >0 


(4.26a) , 


f-^(x)=zk C 


1 /kxtanA e 


/ 2 tt /’THF 


) 


■lyx 




■ix - 3 


-1 


for x >0 

d 5 [l-(l+i) E*(2y5)3 - )E*{ 2 y(-x)}] + 


• (Hi)' 


for x <0 


(4.26b) 


We have been able to calculate f^^(x) analytically for 
n=0,1. Unfortunately, the same cannot be done for the more 
complicated integrals in f— (x), the contributions to £(x) 
from the branch cuts at X: =±K. However, we notice that these 
branch points are weak, and so, their contribution to f( 5) 
for large 5 should be small in comparison to that from 
fi^ ( 5 -*■<»). With this in mind we make the following 

approximation in (4.25b): 

+1/K 


if (X.T) = -1/7 


-yf 1 


(R/(X-C)^ + Y^) 


KY 


-1/K 


/(x-c)H y2 


(5)(4.2T), 


.+ K, 


where we may now replace f— (§) by its small-C behavior. The 

H *★ 

latter may be calculated by deforming C, to C, and C, in 


fig. 4.4 for x>0, x<0, respectively. The integrals which 


I 
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result turn out to be related to the well-known Exponential 
integral, which for small argument (here, Kx=0) has the 


behavior indicated by Bender and Orszag, 

+K 

final expressions for fj^ ( E) are 


(2s) 


pp. 251-252. The 



Vu sin 2^^ 1?* 


n+1 r-. 


n 


.2 


■) [K (-5)]" 




W-5)] - Z Miiui 

j=0 , jl U-n) 

n 


for CO 


^ (K5)\n (K5) -y (-D^r" (Ks)J 

0=0 — j fTT-n) 

05^ 


T 7 ) 


1 

1^" 


(- 1 ) 


n+1 


ni 




i ) (K5)" 


J = 1 


forC>0 (4.28) , 


where Yjj is Euler’s constant, approximately .577. 

Since we are interested in the small-x and large-Y 
behavior of I^(x,Y), we may approximate the Hankel function 


in (4.20) by its form for large argument, 
then becomes 


Equation (4.25a) 


nr 


li'(x,Y) ^ K e"T " 


ikx? , 

dE e 'T (E) 


(4.29) , 


on 


+ K 

and similarly for I- (x,Y), but with limits of integrati 

+K . 

-1/K and +1/K. In (x,Y) we may also make the 

approximation expdKx 5/Y)z1+iKx 5 /Y. 

■ +_y +K 

Performing the integrations in and , and 

substituting the resulting expressions into (4.18), we obtain 
the first two terms in the asymptotic expansion for the 


t 
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I 


I 


I 


I 


I 


loading as we approach the leading edge 

1 1 


- iKY 

n = D + (i-T) §_ 

‘^*2-° z-,A^ / k, - /T v1« 




(^41^ /ti)? _L f (c« + lx c°) 

1” V jiY 

/lscuJnd^ (ii^\ / kxtanA X, 

L 2^ ? ^ vn=n \ 


/T 


Ti^\yH|x 


- 2 +C jx) 

IT Y 


•3 


/=rrior I rnnorr/"^ 

yi^j 5 T 


+ g ^X soun d 


where stands for the acoustic wavelength and 


(4.30) ; 


0 

Ci= - 

' (f n; * 

1 - 

I 

j = 

iil 
= 1 j'J 
00 

^+1) 

(4, 

.31a) , 

Q 

C2 = 

1 (f 


i 

- jl 

ili 1 
7^ (3^ 

(4, 

.31b), 




j 

= 1 

CO 




II 

o 

1(^ J- 

- 

1 - 

I 

j = 2 

1 

j! (j’l) T3+1) 

(4, 

.310, 

1 

Cz = 

i (x+ I - 

3 

|) - 

1 

2 

CM 

•r-y 

1 

1 

TTU-DTTPJT 

(4, 

.31d). 


Figure 4.5a is a plot vs. x/y of the magnitude of the 
first terra in square brackets for n=0 in the asymptotic 
series (4.30) for the correction to the two-dimensional 
loading near the leading edge; the values of M and A are 0.8 
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and 62. Since we interested here in the case of subsonic 
trace speed of the disturbance, in the plot the quantity ^ 
has been replaced by -i*^-U^, as previously discussed. Note 
that the first corrective term is larger on the plate side 
than upstream, as we should expect. Figure 4.5b is a similar 
graph for the second term in square brackets appearing in 
(4,30), corresponding to the n=1 term in the asymptotic. As 
indicated, this terra is to be multiplied by the ratio of the 
acoustic wavelength to the distance from the side edge. 

Next we investigate the behavior of (4.17) near the side 
edge Y=0+. Expanding the Fresnel integral for small 
arguments we have that as Y approaches 0+ 

P*(x,Y>0,0) = (1-1) oqWq. U 

2ttv4^ 



/2 y (-1)’^ 2 nl _ Y 
Vir ^ (2n+iT! 

n=o 




/TPP- 


in+1/2 

(4.32) 


We now call 




(4.33) , 


and recall (4.22) in order to study the contributions from 
the points Ai=+u. It follows that 


+U 

' / 2 kytanAN n-^yz 

i /TiF j 


-iX’iX •• / 

dAj.e (Xi-].i)^(X:tm)^''"^z 2 (4^34). 
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which indicates that lj^^(x>0)=0. For x<0 we deform C, to C?"* 
in fig. 4.3 and find that 


l-\x>0) = 




^2 kytanA^ n+i /2 d(-x}n 

2 (-1)’^ s1n CTr(^+^)] r(iT^) 
fz kytan\\ ^^ 

I 






dn 


d(-xr 


(4.35). 


~+K 

In order to calculate (x) for x>0 and x<0, we deform 
Cl to Cl* and Ci** in figure 4.4, respectively. The results 
are 


ii\x) = (lisiSlii)'' ? ^1+ ’) (4.36a), 


-iKx 


for x>Q 


n 

,+K 

I 


n(x) 


(1+i )/ir 


AytanA 

\n-i/ 2 / 

/P5F 

i 

I 


K+v 


nj/ Y 

(n+1/2) jztTP (^«36b). 


Substituting for I^i^(x) and (x) back into (4.32), we 
obtain the behavior of near the side edge (the n=0 term): 


~+K 


P* = (i-D Chvh U 


,7= 


7HF 


/kx 17, AvtanA 

C(l-i) e'>“- 



0-1) /i^T7 


K-u 


(l-H)/~ ( 
v K+u 

1 


< TT xP / 2 for x<0 

K+u ^ ' 


■sjr for x>0 


(4.37), 
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which has the required behavior for Y=0+ for x>0. For both 
x<0 and x>0 we see that the solution behaves as 1/| xl , where 
X, recall, is normalized by the acoustic wavelength. 

In summary, in this chapter we have obtained an 
approximate solution for the pressure field surrounding a 
quarter-infinite plate passing at high speed through an 
oblique gust. The problem was motivated by the need for a 
better model for the strength of acoustic dipoles near the. 
tip of a blade passing over a vortex when the interaction 
angle is such that the resulting disturbance travels 
subsonically through the still fluid; that developed in- 
chapter 3 overestimated the strength of these tip acoustic 
dipoles. 

We arrived at an expression for the loading on the plate 
(eqn. (4.17)) which we later went on to show satisfied the 
required edge conditions. Also, we showed that away from the 
side edge this result approached asymptotically that obtained 
earlier in chapter 2 for the two-dimensional loading, which 
vanishes upstream and has the 1/yT behavior on the plate; it 
turned out that this was by far the dominant term in the 
calculated asymptotic series. We conclude then that the 
result in (.4.15), from which the loading was calculated by 
setting Z=0, is a good approximation for the pressure field 
for high-reduced frequencies. 

Normally, the next step at this point would have been to 
translate the solution from the present reference frame fixed 
on the plate to one fixed on the ground or still fluid. 
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Unfortunately, the small amount by which we violate the 
upstream boundary condition of P(x<0 , y>0 ,0 ) =0 then becomes 
unfairly magnified by the forward enhancement effect of 
traveling sources, essentially the same (l+Mx^/rg) factor 
present in eqn (3.1^) for that other ground-frame solution. 
The only course left for us now seems to be to extract 
rationally from (4.17) an approximate representation for the 
loading which keeps the physical features we have worked so 
hard to get - the edge behaviors, but that at the same time 
vanishes upstream. This task we undertake in the first part 
of chapter 5 and obtain a relatively simple expression which 
allows us to reintroduce the trailing edge into the model. 
Since so much of the analysis in chapter 5 depends on the 
loading expression in (4.17), we rederive the latter next by 
an independent method, a method based on the idea of source 
cancellations and first used by Landahl to solve unsteady 
transonic flow problems. 


Appendix: Calculation of loading by source cancellations 

(s) 

We seek a load function r (x,Y,0) which when added to 

C) 

p* (x,Y,0), the two-dimensional loading from the leading-edge 

(M (s) 

problem given in eqn. (2.15), we will have p; 4p* 'zO for Y<0. 

The superscript s stands for "side-edge”. The boundary-value 

(s) 

problem for p* is 


XX YY 




Pizz 


= 0 


(4.38) , 
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with 


pi^^(x,y<0,0) = - pP^(x.Y.O) 

pI|^(x,Y> 0 , 0 ) = 0 . 


(4.39a) , 
(4.39b), 


We define the transform pair 


p1^^(Ai,Y,0) = 
P*^^(x,Y, 0 ) = 


|S*^(x,Y,0) 

(4.40a) , 



pI^)(A:,Y,0) 

(4.40b) , 

_<a 

and (4.39a,b). The 

resulting 


^s) 

boundary- value problem for pi is 


p*!! p^!! + -K) = 0 


YY zz 


(4.41 ), 


with 


pi^^ (Ai,Y< 0 , 0 ) = - pP^Ui,Y, 0 ) 
pi^^Xi,Y>0,0) = 0 


(4.42a), 

(4.42b). 


e“) -(s) 

Landahl -has shown that p* (Ai^Y>0) may be obtained by the 


integral 

:(s) 


pi^nAi.Y>0,0) = -1 

IT 




•' 0 


A n^ „-i/K^(-n+Y) 

an P* (Ai,-n,0) e (4.43). 


n-*- 
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(i) 

Substituting for p* , which from eqn. (2.14) we see is (with 
the exp C-ik ^tartV Y/v'l terra included) 


pPUai,Y, 0+) = -OnW. U 


1 


1 




.- IkxtanA Y 

/I 


into (4.43) we obtain that for Y>0 , 


Ui,Y>0,0) = CQV^ _J_ 


/2tt /MF /F 


•{1 - (1-0 E [Y - /K^]} 

/PF 


^ “i kytanA Y 

/iPF 


(4.44) , 


so that for I>0 




U e-lMY^ 
APIT 


/ 1-r 


E[Y ( kxtanA 

i 

(4.45) , 


/1-M 


from which we finally find that 
p 1” + pi=> = f i-l) a,wh U 0-^'^l^t.nAY 

y^F/ k T~'' 5^ 


(4.46) . 


»CD . - 

• dAi _ e' E CY ( kytanA _ )] 

Ai- ir /IPF 

. 

«• 00 
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Chapter 5: Aeroacoustic model for a blade of semi-infinite 
span with improved side-edge loading 

In chapter 2 we saw that the effect of the passage of a 
two-dimensional airfoil through a gust is essentially one of 
diffraction. That is, the pressure field and loading for the 
interaction with a' single gust of wavenumber oib/U is given by 
a continuum of wavenumber solutions - recall eqn. (2.14) for 
the leading-edge effect, for example. 

For the three-dimensional wing with side edge, a similar 
diffraction of the gust spanwise wavenumber k^tanA must occur 
in the y direction. With Z=0+, the Xa integral in eqn. 
(4.16) gives the side-edge diffraction loading predicted in 
our model. We observe that besides the pole at 
\z =kj^tanA/»^ 1 -M^, the integrand also has a branch point at 
Aa . For large Y, the pole term gives the predicted 
spanwise loading the same simple sinusoidal 
exp(-ik^tanAY//1 -M ^) behavior as that of the infinite-span 
wing - the spanwise variation of the input gust upwash; it is 
the only term present in the Y transform of the cut-off 
loading model later developed in this chapter (see eqn. 
(5.21) below). The branch point singularity at Xa - Xi , 
therefore, is responsible for the departure of the loading in 
(4.16) (for Z=0) from its two-dimensional form for arbitrary 
Y. We may say that this is the part of the Xa integrand 
which provides the spanwise diffraction effect for the 
loading: it allows the loading to behave in ‘'magnitude as /~y 
at the side edge and, as we find from the analysis below, 
causes it to have a phase variation near the side very 
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different from that assumed in chapter 3 and in the cut-off 
model of this chapter. Notice that (4.16) also indicates 

that the predicted leading- and side-edge diffractions are 
not independent (even for Z=0). If they were, the integrand 
would be a product of two functions: one depending on 

alone; and the other on ^ 2 , Rather, we observe a coupling. 

The analysis here begins with (4.17). Rationally, we 
make a high-frequency and large interaction angle A 

approximation which essentially eliminates the coupling 
between the leading- and side-edge effects. This allows us 
to obtain an approximate, but closed-form, lifting-surface 
theory for the passage of a quarter-infinite plate through an 
oblique gust. The resulting aerodynamic model predicts the 
strength and phase of tip dipoles more accurately than that 
of chapter 3. Yet, it is still simple enough to allow us to 
calculate the acoustic field in closed form. 

As mentioned above, in this chapter we also develop a 
cut-off aeroacoustic model analogous to that of chapter 3: we 
look at the acoustic field due to a loading pattern given by 
the two-dimensional loading from chapter 2 for y>0, and by 0 
for y<0. 

Approximate lifting-surface theory 

In chapter 4 we developed an. approximate expression for 
the loading on a quarter-infinite plate passing through a 

i • 

gust. This result was derived first by the Wiener-Hopf 
technique and given in. eqn. (4.17), and then again by the 
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method of source cancellations and presented in (4.46). For 
convenience we write it once more. 


/-rfsir + u J”^ /RF 




The integral in this solution may be re-expressed by 
means of the convolution theorem as 


Ci-1) 


dF 

I± g(x-c) 


(5.2) , 


where 


9^^) = £[L. E [Y ( 


' kxtanA 


- A^)] 


(5.3) . 


We now split the integral in (5.3) into three segments, 

as indicated by (5.4); and then for each segment interchange 

orders of integration with the Fresnel integrals. The 

resulting set of integrals is given in (5.5). 

-K K » 


/n g(x) = 


(5.4) . 


Y kytanA 


9(x) = f /TlfF 


Y kxtanA 
rkytan 


/2:rt 


+ /K^- (t kxtanA 
dXe-*^^ 


I Ik^- (_t - kytanA )^ + 

/TW 
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Y kxtan/V 

/TW 


dt e 


it 




kxtanA + 

/TFF 

Y kxtanA 
/Tm^ 

Y kytanA 


lt_e''^ f 


y^TTt 


dxe'<^>^ 


''K^- (t - kytanA )^ 

^ /TW 


-K 




dte 


it 


'0 /Zirt 

Y + i 


dXe 


-iXx 


’ AIF 

« 

Y kxtanA 


dt e 


it 


^t kxtanA^ ? 


■'^TTt 


dXe 


-iXx 


(5.5). 


/T^ 

These may be further manipulated into 


a new set of 

terms, (5.6), where we number each term (1) through (5) for 
ease of reference in the discussion which follows; 


(1) 

2Tr<S(x) E 

Y (kxtanA 
^ /1-Mz •' 

C2) 

-2 sin Kx 

X 

{■ 

Y kxtanA 

I- /1-M^ 

C3) 

2TTd(x) ^ 

' kytanA , 
r ■ i»— — ■ + ICO 


Y ( kytanA 

^ /M=F 


- K 


' dt e'‘ 


Y kxtanA /2Trt 

/TfF 


(4) Y kxtanA 


/T^ 


dt e 


V/ kxtanA >^Tit 

‘V ' , 


it (t . kxtanA )" 


/‘iTit ^ 

-K) _v 

Y /UfF 


die' 
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( 5 ) 


Y kytanA ^ 




dt e 


It 


Y kytanA 


i^ift 


1(2_ (— _ kx tanA ^ 
N Y /TImT 


- l^rifTEZIEK) 


dX e 

T 


-Ux 


(5.6). 


I 


I 


Since £(«)=( 1+i)/2 , we observe that for large Y terras 
(2)-(5) vanish; we are left with (1) only. 

For Y=0+, integrals (3) and (5) may be deformed in the 
complex t plane from the positive real axis to the positive 
imaginary axis (and the contribution from the arc at infinity 
is zero). We then obtain, for terms (3) and (5) respectively 




gi u/4 




dge~^ 
° /2i|~ 



•’o TZttC 


fP+cTF’ 

dxe' 


•lAx 


. 2 tt 6 (x) 


(5.7a), 


(5.7b), 


For Y-j*0+, the inside integral in (5.7a) becomes (x); 
so that the first-order contribution of (3) + (5) to (5.6) 

near the side edge is zero. Higher-order contributions 
vanish as well. 

In order to obtain the behavior of integral (4) for 
Y-J-0+, we integrate by parts as follows: 

t = Y kxtarA 

/T-M^ 

t = Y (iS2LtiOA _K) 



(5.8a), 
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/I /kvtanAY " 2 sinj<^ 

Vire J ^ (5.8b), 


which may be interpreted as the sraall-Y behavior of 

Y (5.9). 

AIF 

We conclude that near the side edge integral (4) cancels 
the first term of (2); so that the behavior of the sum of all 
terras in (5.6) for Y -*■ 0+ is 


iim (4) = 2sin Kx r 
Y-H)+ X 


g(x)s {2tt6(x) + 2 s1n kx} E 

X 


Y^ kx tanA 
/I -M* 



(5.10a), 


from which 


X- C 



• « 

Y^jcxtanA 




(5.10b) . 


We note now .that near the leading edge (x-*- 0), the 
second terra in curly brackets in (5.10b) behaves essentially 
as 2K forE=0. In the convolution (5.2), it contributes a 
terra proportional to the constant K//fuf , which besides being 
small for large interaction angles A (jy] >>K), does not 
contain the 1/ •/T' singularity for .the loading near x=0+, and 
so may be neglected in favor of the 2 iT(S(x-E) term in (5.10b). 

I* ' 

In summary, we find that for both near the side edge, 
and away from it, we may make the approximation 
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I 


I 


I 


i 


g(x-?) = 2ir6(x-C) E 


Y 


ky tanA 


/T^ 


•K) 


Evaluating the convolution (5.2), we finally obtain the 
approximate lifting-surface result for a quarter-infinite 
plate passing through a gust a large interaction angle A 
(subsonic trace speed of disturbances); 


p,(x,y>0.0) = 1 Jl £flWQji e 1 E [y (fcx tanA 




7 




l-M' 


+ u 


^ E [Y (^ 
A 


K)] 


(5.11a), 


which, since E(a = 0) -1^ A indicates that 


p*(x,Y^'‘‘,0)s i 2 pnWn U 

/TFF > 


IT 




,-iyx 

1 J_ 


-iKY 
e /T 


kx tanA -K 


1/ 2 


/1-M 


(5.11b) . 


We find that the Y variation of the loading phase near 
the side is given by exp(-iKY), rather than by 
exp(-ikj^ tanAY'/ /1-M^ 1, as in chapter 3 and in the cut-off 
model developed later in the present chapter. Since the rate 
of decay of the dipole strength near the side is given by the 

iji 

constant (kxtanA-K) , which is >>K for. large A , we notice that 
the spanwise variation of dipole phase is small throughout 
the region we call the "tip region” (fig. 5.1). For large A, 
therefore, the phase of the spanwise loading may be regarded 
as constant for Y=0+. 

The uncoupling of the load into the product of a 
function of x and one of Y allows us now to cancel the 


I 
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solution in the wake by the procedure used in chapter 2 for 
the trailing-edge problem. As we would expect at this point, 
the expression for the pressure distribution on the 
semi-infinite span blade with' side edge turns out to be 



which becomes the two-dimensional loading far from the side 
I edge as it should. 

Acoustic field for a semi-infinite span blade with side edge 
passing through an oblique gust 

' As we noted above, the simplified loading expression 

allows us to solve a new trailing edge problem; the details 
of this analysis are no different from those in chapter 2 and 
so will be omitted here. We now define a pressure 
P(x , Y, Z, t)=;^(x, Y, Z, t)exp( iut+ik^ M2x/( 1-M2 ) ) , which is a 
solution of 

\ 

j 

^xx * P*YY * P*ZZ * 

1 

subject to 


p*(x,Y,0) = p* (x,Y,0) 


(5.14) , 
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the load we have calculated above for a semi-infinite blade 
passing through a gust. The solution procedure of this 
boundary-value problem is similar to that in chapter- 3 for 
the rectangular rotating blade. Only one new integral comes 
up which we did not have before: the Fourier transform of a 
Fresnel integral; 


dY Y ^ ry tan A -K) 




/PF 


(5.15) , 


which exists , for ImCVj ) >Im(kj^ tanA/Zl-K'^ ) and may be readily 
calculated either by integration by parts, or by 
interchanging orders of integration. Thus, we obtain that it 


is equal to 


3Tn’ 

^ g 4 AxtanA -K 

J /TiF 

L kytanA\ /r ;r 


( 5 . 16 ) , 


At this point we allow r^= /x^ +y , and obtain the 
far-field pressure by the method of stationary phase. The 
procedure is completely analogous to that presented in detail 
in chapter 3. The final result for the acoustic signal 
P(n,y,(fi, t) perceived by an observer on the ground is 


P(r3,Y>4>>t) = -i pnWn U M 
27 t : 


_M lw_(t-bra/Ca) 

.. (y.cj,) 

y/j^ >"a ( l+Msln(j)COSY) 

X 


(5.17) , 
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where the ”/y'» subscript on D/-(y» 4>) stands for /T^load 

''y 

behavior at the side edge. Dy-(y,4)) was found to be 


D _ (y,4i) = 
Vy 



the three-dimensional function which describes the acoustic 
directivity pattern in terms of the spherical angles y , in 
fig. 5.2. 


Notice that the blade's tip is located at the origin, 
so that its semi-infinite span extends out from Y=0. This is 
in contrast to the system we had in chapter 3, where the 
blade tip was rightfully "outboard." Still, here as in 
chapter 3 the interaction angle A is such that the trace 
velocity of the unsteady blade loading is in the positive Y 


direction 



Acoustic field for a semi-infinite blade with square cut-off 
loading 

Here we write down the solution of the boundary-value 
problem for p^(x,Y,Z): 


with 


P* ■*' 
XX 


p* 


YY 


p* 


+ K2 


ZZ 


P* 


= 0 


(5.19), 


^*(x,Y.O) = P* 2 _d (x.Y.O) forY>0 (5.20a), 
p* (xJ,0) = 0 for Y<0 (5.20b); 

where p* (x,Y,0) is given in eqn. (3.4), 

2-0 

In the solution process for the present problem, we find 
that the spanwise Y transform of the loading reduces to the 
evaluation of the simple integral 


dY ) Y ^ 


i 


CA2- kxtanA ) 


(5.21), 


which converges for Im(A 2 ) >Ira(k^tanA/V ). Multiplied by 
(1+i)/2, this result for the cut-off solution replaces that 
in ( 5 . 16 ) for the spanwise spectrum of the blade loading with 
improved side-edge behavior. 

The acoustic far field for a semi-inf indt'e span wing 


with square cut-off. of the loading at the side edge is 
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ifiLlli&L 

P MO-MiW-* 1 e l+Msin<DcosT 


P -^ ^nWo_U MiMllV** _L ■ Dc.o(r,^) .5 22) 

2^^/TFF / kv"'7^ (l+H+in4.cos )" * ^ 

V 1“M 

where the "c.o." subscript in Dc,o^.Y>'l>^ denotes "cut off." 
was found to be 




( 5 . 23 ) . • 


Predicted directivity patterns for D and D 

^ c.o. 

The examples we show next of D ^ and ^ focus on the 

acoustic effect of satisfying the side-edge loading 

1 

condition. We do not, therefore, investigate here again the 

effect of the trailing edge on radiation (see fig. 3.8), 

' because this is a chordwise effect. Nor do we look for 

, $ ■ 

changes due to variations of Mach number - although most 
compressibility effects cannot be extricated from those of 

\ 
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frequency, we saw in chapter 3 that changes only in the Mach 
number affect the forward enhancement of the signal; so that 
we may choose to regard M as really a chordwise, or 
flight-direction parameter. 

Here we investigate the acoustic effect of reduced 
frequency k , which affects the phase and magnitude of the 

A 

loading near the side; and of A , which besides also affecting 
the side-edge loading, controls the trace speed U/tanA of 
disturbances along the span and the absolute speed lI/sinA of 
load fronts through the still fluid for a given M. If we 
change A and keep M constant, we may compare the degree of 
cancellation of signals radiated by the semi-infinite blade 
span for the blade with improved side-edge loading, to that 
for the blade with cut-off loading. 


So that we may also compare D fjj) to those for the 

rotating wing of chapter 3, we have plotted the quantity 


D* / tanA- Ji_ (y.4>) 

y5T Z/iT' v /] (1+H sinfjjcosY)^ 

in figs. 5. 3-5. 5; and for a similar, reason, the quantity 

D*co.CY.<f) = 

in figs. 5. 6-5. 8. As in chapter 3, every figure shows the . 
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above functions in five planes, each corresponding to a fixed 
value of the azimuth ? • 

From figs. 5.3 and 5.6 for the correct side-edge loaded 

blade and the blade with cut-off loading, respectively, we 

note that I? is generally greater than on the planes 

shown in parts (a)-(c) of the figures. D (y =1 50®, 0<(|)<ir/2) in 

fig. 5.3e, however, takes on values larger than those for 

D* (y =1 50®, 0<<i»<7r/2) in fig. 5.6e. This perhaps unexpected 
c.o. 

peak in D*^y =1 5Cf,4> =79“- ) can be attributed to a local 
reinf orceraent of acoustic signals - recall that the spanwise 
variation of the phase of tip dipoles is different for the 

models from which D* and „were calculated. 

c.o. 

As we noted earlier, we would expect the differences in 
the directivity patterns between the two models to increase 
as A increases, because even though the region near the tip 
where dipole strengths differ substantially is reduced in 
extent, the differences in spanwise variation of the phase of 
tip dipoles are then larger. This is in fact the case, as we 
may see from figs. 5.4 and 5.7: the differences in radiation 
patterns between cut-off and improved side-edge loadings are 
greater for A =80® than for A=62*J Hotice also that for both 
types of loading the overall radiation for a= 8CP is 
substantially reduced. As in chapter 3, we attribute this 
reduction to the fact that the effective Mach number of load 
fronts M/sinA far from the side edge is lower for a=80® than 

O ' 

for A=62; so that in both loading cases there is then more 
cancellation in the acoustic field radiated by the 


! 
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serai- infinite span of the blade. We conclude that the 

differences due to changes in M/sinA in spanwise loading 

phase between the correct and cut-off models may cause 

significant differences in the predicted acoustic fields. 

By varying A while keeping ub/U and M constant, we 

* * 

observed differences between and D essentially due to 

the different way each of the two acoustic models reacted to 

a change in the absolute Mach number M/sinA of load fronts 

far from the side edge. But the tip-loading expressions 

derived earlier indicate that a change in K (by changing 

ojb/U, with M/sinA constant), also affects the differences in 

spanwise variation of tip-loading phase between the two 

models, and should, therefore, also affect their acoustic 

* 

fields. Figs. • 5.5 and 5.8 show and , respectively, 
for ojb/UslO while keeping the same M and A values of figs.' 
5.4 and 5.6. They indicate that the differences between 
and Dc.o.fof* ub/U = 10 are not much greater than those for 

ub/U=4; so that pure frequency effects on spanwise phase 

differences between the two models have a lesser acoustic 
effect than does a change in M/sinA. 

Finally, we observe from (5.17). and (5.22) that the 
acoustic models developed in this chapter share with that of 
chapter 3 the 1/k behavior for large k . This means that 

X X 

they all predict acoustic signals which are qualitatively 
similar for small times. We note, however, that for ^ =0 

(large times) the models developed in this chapter for a wing 

of semi-infinite span are not well behaved. The physical 


I 
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reason, for the large-time breakdown of the theories is that 
an infinitely long time is required for pulses radiated near 
y=oo to reach the observer, and so, to contribute to the 
signature which he hears. This issue of convergence, in time 
did not arise in chapter 3 because the blade there was of 


finite span 
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Chapter 6: Summary and Conclusions 

We have developed here three basic aeroacoustic. models 
f-or noncompact blade-gust interaction in subsonic flow. We 
now briefly summarize our results and draw some conclusions. 

Infinite-span airfoil passing through a gust 

In chapter 2 we found that frequency and compressibility 
effects combine into a single parameter which sets the 
criterion for chordwise noncorapactness . We observed that an 
airfoil responds to an input single-wavelength gust 
essentially as an aerodynamically diffracting surface - the 
solution has a continuous wavenumber spectrum. The 
closed-form expression we have derived for the acoustic field 
in a laboratory frame indicates that tone radiation patterns 
from compact, and noncompact airfoil loadings may differ 
dramatically. The complicated lobing present in the 
noncompact case is caused by the reinforcement and 
cancellation of signals from a chordwise distribution of 
dipoles with large differences in phase. The unsteady 
loading for an airfoil which is acoustically compact may be 
represented by a simple dipole, which in turn radiates a 
simple dipole acoustic field. 

We verified the expected lack of acoustic field for 
interaction angle-Mach number combinations which give load 
fronts a propagation speed less than that of sound. 



Rotating blade of finite span passing through a gust 

In chapter 3 we used chapter 2’s two-dimensional 

solution in a superposition scheme to model radiation from a 

rotating blade of finite span. We found that the acoustic 

field could be calculated only for cases with subsonic 

disturbances; the impossibility of extending the analysis to 

include supersonic trace speeds is an artifact of the 

(io) 

solution procedure, which is based on Landahl's edge 
independence concept. We studied Mach number, frequency, and 
interaction-angle effects on the radiation pattern of a 
single tone. The results were presented on the ground, or in 
the acoustic frame of reference. We refer the reader to the 
discussion at the end of that chapter for detailed 
conclusions. We found that for a given set of physical 
parameters describing the wing-gust interaction, ie. M,A , 
and gust amplitude w^, the noncorapact theory predicts that 
the tone amplitude at every field point decays algebraically 
as where is the reduced frequency. This is in 

contrast to the high-frequency behavior of the tone amplitude 
predicted by compact theories (se.e, e.g. ref. 3); these 
predict that for given values of M, A, w at high frequencies 
the tone amplitude ’’levels off” asymptotically to a constant 
value. 

We applied the leading-edge part of the theory to 

» • 

predict the transient signal emmitted by a finite-span blade 
passing over a potential vortex. We found that the level of 



the signature was a sensitive function of spatial position. 
We attributed this sensitivity of the total signal to 
cancellation and reinforcement effects of the tones making it 
up - another consequence of having a noncompact loading. 

Semi-infinite rectangular blade with improved side-edge 
loading 

Although we realized in chapter 3 that the model 
developed there overestimated the strength of dipoles near 
the blade tip, we had no way of knowing how badly we- were 
modelling the spanwise variation of the loading phase. In 
chapter 4 we calculate an approximate solution for the 
three-dimensional pressure field surrounding a corner plate 
passing through a gust. We find that for large interaction 
angles the spanwise component of. tip-loading phase is nearly 
constant. The actual variation of the tip-loading phase 
implies a different pattern of reinforcement and cancellation 
for signals in the far field. 

In chapter 5 we simplify further the aerodynamic result 
of chapter 4; we then use the new result to calculate the 
acoustic field for a semi-infinite span rectangular blade, 
with improved side-edge loading, passing through a gust. We 
also calculate the acoustic field for a similar blade, but 
with loading given by the two-dimensional response of chapter 
2 abruptly cut off at the tip. This cut off loading is 
analogous to that used in chapter 3 for a finite blade in 
that the spanwise loading contains only the spanwise 
wavenumber of the input gust. We find that the directivity 
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patterns differ mostly when A , the interaction angle, is 
large. The difference can be explained by noting that the 
two models then differ greatly in the -way they predict the 
phase of tip dipoles. 

We observed that the three-dimensional models developed 
in chapter 5 share with that of chapter 3 the 1/k^ decay of 
the tone amplitude as a function of high-reduced frequency. 
As discussed above in chapter 3's summary, this algebraic 
falloff indicates a fundamental difference in the way 
radiation patterns are predicted by noncompact and compact 
models; for given values of M, a » and Wq the latter predict 
that at high-reduced frequencies the amplitude of a tone 
reaches a constant value. 

Conclusions 

The signal from a blade which radiates as a noncompact 
acoustic source probably has a maximum peak value which is 
smaller than that for a similar signal from a blade which 
radiates as a compact source. This conclusion is based on 
comparison of the high-frequency behavior of the amplitude of 
tones which make up the acoustic signal as predicted by 

n 

compact theories developed elsewhere with that predicted by 
the noncompact theories developed here. Unfortunately, 
because the predicted signals in ref. 4 were calculated for 
supersonic trace speeds only, direct comparison of those 
results to those computed here is not possible. 

As noted above in the summaries of chapter 3 and 5, for 



given values of M, Wq , we expect a rapid decrease of the 
level of high-frequency components for a noncompact 
situation. For given values of M, a, compact theories 
predict a spectrum which reaches a constant value at 
high-reduced frequencies, thus suggesting that tone 


reinforcement could result in a louder signal. 
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Fig. 1.1 (a) Relative positions of blades artd vortices for a 
single-rotor helicopter; (b) Stages of a blade-vortex 
interaction; (c) unsteady blade force and its time 
derivative. 
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Fig. 2.1 The passage of an infinite-span wing through an 
oblique gust. 
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Fig. 2.2 Regions of analy ticl ty O i © complex contour 




t 


98 



Fig. 2.3 Deformed contour C*, for x>0. 
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Fig. 2.4 Deformed contour C**, for x<2. 
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Fig. 3.4a-e Plot of the mag 
acoustic directivity 
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Fig. 3.5a-e Plot of the raagnituae of the three-dimensional 
acoustic directivity for k^ = 10j M=.8,A =62P, L/b = 15. 
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Fig. 3.10 Part of the directivity pattern of the peak signal 
as predicted by eqn. (3.29). ; " 
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Fig. 3.12 Predicted spectra of signatures at positions 
indicated for h/b=.5 and h/b=.1. 
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Fig. 4.2a,b Regions of 0 and 0 analyticity in the complex 
ai,A 2 planes. Integration contours Ci, Cj are indicated. 
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Fig. 4.5a,b Plots of first, and second terms in brackets in 
.eqn, (4,30). 
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Fig. 5.la,b Spanwise loading predicted by (5.11a,b): 
magnitude and phase for Aa i* ground reference 

frame. 
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Fig. 5.7a-e Plot of the magnitude of the' three-dimensional 

acoustic directivity D* ( Y, <J>) for k =4 , M=.8, ifc80°. 

c.o . 
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Fig. 5.8a-e Plot of the magnitude of the -• three-dimensional 
acoustic directivity for k^ = 10, M=,8, a=62®. 
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